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CHAPTER  I 


INTRODUCTION 


This  report  represents  a  compilation  of  the  important  as¬ 
pects  of  an  analytical  and  experimental  study  concerned  with 
non-uniform  plasmas.  The  information  presented  in  this  report 
will  not  necessarily  represent  the  chronological  development  of 
the  experimentation,  but  rather  a  summary  of  the  physical  prob¬ 
lem  and  the  methods  of  solutions. 

The  experimental  work  was  done  using  a  shock  tube  to  pro¬ 
duce  the  plasma.  A  detailed  description  of  the  shock  tube  facil¬ 
ity  is  presented  in  Chapter  II.  Of  particular  interest  is  a 
description  of  a  mechanical  valve  which  has  been  successfully 
used  to  replace  the  conventional  driver  and  diaphragm  sections 
of  the  shock  tube.  In  addition  a  description  is  given  of  a  thin 
film  resistance  gauge  and  amplifier  circuit  which  has  been  found 
to  be  a  useful  trigger  under  a  wide  variety  of  conditions. 

An  integral  part  of  the  work  performed  on  this  contract  has 
been  the  development  of  a  spatially  resolving  electron  density 
microwave  diagnostic  technique  useful  in  a  range  of  electron 
densities  infrequently  explored  by  other  investigators.  This 
technique  permits  accurate  spatial  resolution  in  plasma  gradients 
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for  electron  densities,  orders  -  f  magnitude  below  that  possible, 
usinq  a  conventional  microwave  interferometer.  After  an  initial 
calibration  phase  in  which  argon  was  the  driven  gas  in  the  shock 
tube,  this  microwave  technique  has  been  used  to  measure  electron 
densities  and  electron  density  profiles  behind  normal  air  shocks 
in  the  range  of  Mach  number  8.5-11  and  initial  pressure  condi¬ 
tions  of  —  1  mm  Hg. 

The  microwave  instrumentation  for  this  technique  is  de¬ 
scribed  in  detail  in  Chapter  III.  The  analytical  formulation 
correlating  the  measured  properties  to  ascertain  the  electron 
density  is  shown  in  Chapter  IV.  Included  in  this  discussion 
are  the  details  of  a  data  processing  program  used  to  convert 
the  raw  data  taken  from  the  microwave  instrumentation  into  a 
usable  form  for  computation. 

When  the  results  of  the  measurements  were  compared  with 
the  values  of  electron  densities  predicted  on  the  basis  of  the 
accepted  semi-analytical  rare  chemistry  program,  an  order  of 
magnitude  discrepancy  was  found.  Specifically,  the  measure¬ 
ments  indicate  an  electron  density  behind  the  shock  which  tends 
to  the  equilibrium  value  at  a  rate  which  is  about  an  order  of 
magnitude  faster  than  the  one  predicted  by  the  numerical  calcu¬ 
lations.  Since  the  shock  conditions  where  the  measured 
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discrepancy  exists  are  relevant  to  the  150,000  ft  altitude  range 
aerodynamic  calculations,  it  was  deemed  pertinent  to  explore  the 
reasons  for  the  observed  discrepancy.  To  this  effect  an  exten¬ 
sive  program  of  measurements  was  initiated  with  different  gas 
compositions  and  the  sources  of  errors  in  the  measuring  technique. 

Chapter  V  presents  a  summary  of  the  theoretical  work  con¬ 
ducted  in  conjunction  with  the  experimental  program.  A  brief 
description  of  the  generalized  properties  of  non-uniform  plasmas 
is  given  and  referenced  to  the  detailed  formulations.  A  section 
is  presented  where  the  diffusion  effects  in  a  non-uniform  plasma 
are  outlined  quantitatively.  On  the  basis  of  the  diffusion  cal¬ 
culation,  the  electromagnetic  properties  of  the  plasma  contained 
in  the  shock  tube  where  radial  non-uniformity  exists  are  dis¬ 


cussed  . 
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CHAPTER  II 
THE  FACILITY 


A-l  INTRODUCTION 

The  explicit  purpose  of  the  following  detailed  description 
of  the  GASL  combined  shock  tube  and  waveguide  facility  is  to 
enable  other  investigators  to  duplicate  any  or  all  parts  of 
this  highly  successful  experimental  arrangement.  Since  it  is 
the  details  rather  than  the  basic  principles  that  are  time  con¬ 
suming  in  the  fabrication  of  such  a  facility,  great  effort  will 
be  made  to  present  all  salient  features. 
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A“2  BASIC  SHOCK  TUBE 

In  this  experimental  facility,  the  driven  section  of  the 
shock  tube  is  also  used  as  a  circular  waveguide.  Because  of 
this  duplicity,  some  precautions  must  be  taken  as  to  the  choice 
of  material,  the  surface  finish,  and  the  method  of  connecting 
tubing  sections. 

A  compromise  material  for  the  tube  is  seamless  stainless 
steel  pipe.  W  \le  this  material  is  certainly  not  the  best 
electrical  conductor  as  would  be  desirable  from  the  waveguide 
point  of  view,  it  offers  the  c  "  'ntages  of  strength  machinabil- 
ity,  corrosion  resistance  and  availability.  The  stock  material 
used  for  the  driven  sections  of  the  shock  tube  were  commercially 
available,  extra  heavy,  three- inch  stainless  steel  pipe.  This 
tube  has  a  nominal  inside  diameter  of  2.9"  and  outside  diameter 
of  3.5" . 

Highly  desirable  from  the  waveguide  aspect  is  that  the 
tube's  inner  diameter  be  constant.  For  use  as  a  shock  tube, 
the  inner  surface  should  be  smooth  with  no  wall  discontinuities 
at  the  tube  section  junctions.  To  satisfy  these  conditions, 
each  ten  foot  length  of  pipe  was  reamed  and  honed  by  a  special¬ 
ist  in  this  field.  The  resultant  inside  diameter  for  all  four 


driven  sections  was  measured  to  be  2.917"  ±  .001". 
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The  problem  of  coupling  the  tube  sections  is  summarized 
as  follows:  for  use  as  a  ./aveguide,  good  and  complete  electrical 


contact  must  be  maintained  at  each  junction,  the  sections  should 
join  with  a  minimum  angle  between  their  axes,  discontinuities 
in  the  wall  should  be  minimized,  and  the  junction  must  be  pres¬ 
sure  and  vacuum  tight  and  mechanically  surong  Also  desirable 
is  a  connection  that  can  be  easily  removed  and  reassembled, 
since  the  insertion  of  various  test  sections  may  be  required 
at  any  junction.  In  addition,  each  section,  should  be  able  to 
be  rotated  and  be  locked  independently.  This  requirement 
originates  from  the  fact  that  the  microwave  radiation  is  polar¬ 
ized  and  both  parallel  and  perpendicular  measurements  are  some¬ 
times  desirable. 

To  satisfy  these  requirements,  the  coupling  assembly  shown 
in  Fig.  A- 2.1  was  devised.  In  this  scheme,  one  end  of  each  tube 
section  is  machined  with  a  right  hand  thread  and  cylindrical 
surface  concentric  with  the  inner  axis.  On  the  other  end  of 
each  section,  a  left  hand  thread,  an  identical  cylindrical  sur¬ 
face  and  a  partical  "0"  ring  groove  are  machined.  This  "0" 
ring  forms  the  seal  between  sections,  and  is  held  in  place  by 
the  inner  surface  of  the  coupling  nut.  These  are  relatively 
simple  machining  operations,  and  if  done  in  the  order  stated. 
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concentricity  and  diametric  tolerances  are  maintained  without 
the  usual  distortion  problems  associated  with  weld  type  coupling 
structures. 

Bronze  is  used  for  the  coupling  nut  ro  prevent  galling  of 
the  surfaces,  and  is  available  as  standard  unfinished  bronze 
bearing  stock.  If  the  tolerances  of  the  cylindrical  surfaces 
are  held  to  .002",  the  maximum  angular  deviation  between  sec¬ 
tions,  excluding  bending,  is  approximately  .1  degree.  Great 
care  and  patience  should  be  exercised  in  the  assembly  of  tube 
sections  in  order  to  prevent  bending. 

This  coupling  arrangement  satisfies  all  of  the  conditions 
previously  mentioned.  However,  in  order  to  assemble  and  dis¬ 
assemble  the  tubes,  the  downstream  sections  must  be  able  to 
slide  on  their  mounting  pads.  Since  the  reaction  forces  on  the 
shock  tube  exert  only  axial  stresses,  there  is  no  need  to 
rigidly  restrain  the  tube  along  its  length,  so  long  as  the  mam 
reactionary  forces  are  absorbed  at  the  diaphragm  location. 

Since  this  is  not  only  feasible  but  desirable  from  the  point 
of  view  of  simple  electrical  isolation  from  ground,  the  sup¬ 
ports  for  the  driven  section  need  carry  only  the  compressional 
weight  load  of  the  tube.  This  being  the  case,  no  extensive 
foundations  or  structure  are  required.  Since  perfect  alignment 
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of  these  supports  to  the  accuracy  dictated  by  the  coupling  is 
not  practically  attainable,  an  adjustable  set  of  "vee"  jaws  can 
be  mounted  on  simple  "I"  beam  structures  to  support  the  entire 
driven  tube.  Such  a  set  of  jaws  are  shown  in  Fig.  A-2.2. 

These  jaws  allow  independent  adjustment  of  height  and 
horizontal  position.  The  vertical  adjustment  is  achieved  by 
rotating  the  central  screw  which  has  right  and  left  hand  threads 
operating  through  the  respective  threads  in  the  "V"  jaws.  The 
vertical  load  on  the  javrs  is  transmitted  directly  to  the  mount¬ 
ing  plate  on  which  they  slide.  Horizontal  adjustment  is  made  by 
moving  the  angle  bracket,  containing  the  center  slot  guide 
bearing,  in  rhe  range  provided  by  the  slots  in  the  bracket. 

The  jaws  are  finally  clamped  to  the  angle  bracket  which  is  also 
slotted  in  the  horizontal  direction.  Satisfactory  support  is 
accomplished  with  two  sets  of  jaws  for  each  ten  foot  section  of 
tube.  With  this  basic  mounting  and  coupling  arrangement  the 
entire  driven  tube  can  be  assembled  in  one  day  with  the  help  of 
a  precision  level,  an  inexpensive  transit,  suitable  "vee"  blocks 
which  will  fit  on  the  cylindrically  machined  part  of  the  tube 
section,  and  a  little  patience  and  luck. 
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A -3  FAST  OPENING  MECHANICAL  VALVE 

General  Characteristics 

This  valve  has  been  developed  and  constructed  to  replace 
the  diaphragm  in  a  small  shock  tube. 

The  conditions  of  operation  are: 

Driver  pressure:  100  to  500  psi 

Driver  temperature:  up  to  200° C 
Opening  area:  7  sq.  inches 

Seal:  vacuum  tight 

Additional  conditions  which  should  be  satisfied  are  re¬ 
mote  operation  and  minimum  of  adjustment  or  replacement  of 
parts.  The  adopted  solution  has  only  one  moving  part,  a  piston, 
that  provides  two  types  of  seal:  an  "O"  ring  seal  for  vacuum 
and  a  sleeve  valve  for  fast  action. 


reverses . 
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(c)  Vacuum  seal  opening:  When  the  piston  starts  to  move 
towards  the  left  the  force  pushing  it  inci eases  sud¬ 
denly  because  the  pressure  starts  acting  over  the 
surface  between  vacuum  seal  and  sleeve  seal. 

(d)  Fast  opening:  When  the  sleeve  goes  out  of  the  driver 
chamber  the  gas  starts  to  flow  from  driver  to  driven. 
The  piston  has  been  accelerated  along  a  distance 
equal  to  the  length  of  the  sleeve. 

(e)  Stopping:  The  piston  compresses  the  gas  left  in 
chamber  A,  stopping  and  oscillating  around  a  posi¬ 
tion  which  is  a  function  of  the  final  pressure  of 
the  tube. 

(f)  Closing:  After  firing,  the  driver  and  driven  are 
evacuated;  the  gas  left  in  chamber  A  expands  and 
closes  the  valve. 

(g)  Charging:  The  gas  is  fed  first  into  chamber  A, 
closing  the  valve  if  step  (f)  fails  to  do  so,  and 
then  the  gas  goes  into  the  driver  through  a  check 
valve. 

The  two  main  facts  to  be  considered  in  the  action  of  the 


valve  are: 
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(1)  The  valve  behaves  as  an  unstable  system  near  the 
firing  pressure.  The  force  pushing  the  piston  jumps 
from  zero  to  a  maximum  value  in  a  very  short  time 

(<  10"4  sec) .  This  fast  action  cannot  be  obtained 
by  controlling  the  flow  of  gas  with  conventional 
valves . 

(2)  The  fast  seal:  Only  after  the  piston  receives  its 
maximum  velocity  is  the  communication  between  driver 
and  driven  established.  For  a  given  acceleration  of 
the  piston  the  opening  time  goes  down  with  the  square 
root  of  the  length  of  the  sleeve. 

The  fast  seal  is  not  tight  and  the  possible  leakage  be¬ 
tween  the  opening  of  both  seals  establishes  a  compromise  between 
length  of  sleeve,  acceleration  and  quality  of  seal. 

Computation  of  the  Main  Characteristics 

The  calculations  are  very  easy  if  we  disregard  secondary 
effects. 

The  main  parameters  are: 

S  area  of  piston 

A 

Sg  area  of  the  vacuum  tight  seal 

transversal  section  of  driven 
D 

P  ores sure  in  the  chamber  A 

A 
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P  pressure  in  the  driver 

B  ^ 

W  weight  of  the  piston 


e  length  of  the  sleeve 

The  pressure  in  the  driven  is,  in  normal  working  condi¬ 
tions,  on  the  order  of  a  few  millimeters  of  mercury  and  we  dis¬ 
regard  it  in  the  ca3 cu' ations.  When  charged  it  is  P^  =  P^,  and 
the  force  closing  the  seal  is  F  =  PgSg.  The  chamber  A  is 
vented  towards  atmospheric;  P^  remains  constant  and  P  goes 
down  reaching  the  firing  pressure  Pf  corresponding  to  zero 


force  acting  on  the  piston;  Pf  =  P^ 


i .  h 
Sa 


if  P 


to  atmospheric  pressure,  the  corresponding  PB  is  the 


is  equal 
minimum  at 


SB 

which  the  valve  can  operate  for  a  given  ratio  —  . 

When  the  piston  starts  to  move,  the  force  pushing  it  back 

jumps  to  F,  =  P„(S  -SJ  -  P.S.?  because  of  the  volume  of  the 
f  BAD  AA 

chamber  A  we  can  assume  that  during  the  fast  seal  open¬ 

ing  of  the  valve,  so  F,  =  P_(S_-S_).  With  this  force  we  can 
^  f  B  B  D 

compute  the  acceleration,  velocity  and  opening  time. 

We  consider  that  the  valve  is  open  when  the  passage 
area  of  the  gas  through  the  fast  seal  equals  the  transversal 


area  of  the  driver. 

The  corresponding  values  for  the  valve  in  operation  are: 
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=  60  psi 

PB  =  300  psi 

S  =16  in2 
A 

SB  =  13  in3 

SD  =  7  ir'3 

=  1600  pounds 

opening  speed  =  120  ft/sec 
opening  time  =  300  microseconds 
Figure  A-3.2  presents  the  calculated  shock  Mach  numbers  as 
a  function  of  the  pressure  ratio  between  driver  and  driven  for 
air.  This  calculation  assumes  an  infinitely  fast  opening  dia¬ 
phragm  and  no  energy  loss  mechanisms  occurring  in  the  flow. 
Experimental  points  (+)  have  been  superimposed  on  the  plot  indi¬ 
cating  the  behavior  obtained  using  the  mechanical  valve.  The 
agreement  equals  or  exceeds  that  obtainable  using  a  bursting 
diaphragm . 


Shock  Mach  Number-  versus 


Solid  line  dervied 


U  M.V  £  r 
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A- 4  SHOCK  TUBE  PLUMBING 

The  governing  factor  in  the  engineering  of  the  shock  tube 
plumbing  is  safety.  Since  the  driver  gas  normally  used  ir  the 
facility  is  hydrogen,  the  possibility  of  explosion  either  in¬ 
ternally  or  externally  due  to  a  leak  is  omnipresent.  While  the 
basic  tube  is  rated  for  over  1000  psi,  the  various  detectors 
and  fittings  are  in  general  not  designed  for  these  pressures. 
Hence,  even  a  limited  internal  explosion  is  serious. 

The  replacing  of  the  conventional  diaphragm  system  with 
the  mechanical  valve  has  considerably  reduced  the  hazards  of 
operation.  With  the  valve  the  use  of  the  shock  tube  has  become 
a  closed  operation.  Since  the  tube  need  not  be  opened,  the  only 
air  entering  the  system  is  introduced  during  the  controlled  fill¬ 
ing  cycle  of  the  driven.  This  air  charging  can  only  be  per¬ 
formed  if  the  driver  is  not  charged,  the  driven  previously 
evacuated,  and  then  the  air  charge  cannot  exceed  one  atmosphere. 
The  driver  cannot  be  filled  until  the  pressure  in  the  driven  is 
again  reduced  below  50  mm  Hg.  Once  the  driver  has  been  filled 
the  only  operation  that  can  be  performed  is  either  "fire"  or 
"purge" .  After  firing,  a  pressure  switch  in  the  driven  starts 
the  purging  cycle  automatically,  venting  excess  pressure  over 
atmospheric  to  the  atmosphere.  Then  the  vacuum  system  is 
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opened  and  the  tube  evacuated.  While  it  is  possible  to  bypass 
these  automatic  operations  to  obtain  unusual  conditions,  it  is 
generally  during  the  routine  operation  of  such  a  system  that  a 
careless  e^ror  is  serious.  Specific  details  of  the  plumbing 
system  are  not  warranted  hare  due  to  the  fact  that  the  use  of 
the  mechanical  valve  instead  cf  a  diaphragm  has  permitted  addi¬ 
tional  simplifications.  The  basic  system  to  be  installed  is 
outlined  below. 

The  combined  driver-mechanical  valve  as  previously  dis¬ 
cussed  need  have  only  one  solenoid  type  filling  valve,  and  one 
small  mechanical  vacuum  valve  to  be  used  only  if  the  assembly 
is  to  be  opened  for  cleaning  or  servicing.  Another  small 
solenoid  valve  may  be  used  if  the  possibility  of  a  high  pres¬ 
sure  dump  is  desired.  The  driven  tube  is  equipped  with  two 
semi-automatic  valves.  The  first  is  an  over  pressure  valve 
which  will  seal  vacuum  but  will  open  a  two- inch  port  if  the 
static  pressure  in  the  driven  exceeds  one  atmosphere  absolute. 
This  valve  mechanically  interlocks  to  the  vacuum  valve  so  that 
if  the  over-pressure  valve  is  open,  the  vacuum  valve  must  be 
closed.  The  vacuum  valve  is  pneumatically  operated. 

With  this  system  a  typical  operating  cycle  would  be  as 
follows.  The  driven  tube  is  pumped  down  by  opening  the  pneu¬ 


matic  vacuum  valve.  The  tube  is  filled  with  the  desired  test 
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gas  which  is  limited  in  pressure  by  the  t.ver-pressure  valve. 

The  tube  is  now  evacuated  to  the  desired  pressure.  The  driver 
is  then  filled.  It  may  be  noted  that  if  the  driver  had  been 
inadvertently  filled  and  fired  during  any  of  the  above  opera¬ 
tions  all  systems  including  the  vacuum  pump  would  have  been 
protected  by  the  over-pressure  valve.  When  the  firing  sequence 
is  started,  a  relay  is  locked  so  that  the  next  pressure  rise 
occurring  in  the  driven  will  open  the  driven  vacuum  valve. 

After  firing,  if  the  pressure  exceeds  one  atmosphere,  the  opera¬ 
tion  of  the  vacuum  valve  is  blocked  until  a  safe  pressure  in 
the  driven  is  reached,  after  which  the  system  will  be  evacuated 
leaving  it  in  a  safe  condition  and  ready  for  the  next  test.  For 
the  above  operation,  it  was  assumed  that  the  main  mechanical 
valve  operating  in  place  of  the  diaphragm  closed  after  firing. 

In  the  event  that  a  malfunction  occurs,  no  further  operation  of 
the  shock  tube  is  possible  until  the  trouble  is  corrected,  but 
neither  are  any  safety  problems  present. 

Two  basic  static  pressure  gauges  have  been  found  to  be 
adequate  in  the  system.  To  monitor  the  driver  pressure  a 
standard  250  psi  Heiss  gauge  is  used  without  valves  or  restrict¬ 
ing  orifices.  An  effective  and  accurate  gauge  for  measurement 
of  pre-firing  pressure  in  the  driven  is  the  Alphatron  gauge. 

This  gauge  (an  old  model)  has  three  linear  ranges,  0-100  microns. 
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0-1  mm  Hg,  and  0-10  mm  Hg,  with  a  zero  shift  raising  the  upper 
range  to  20  mm  Hg.  Supplementing  the  Alphatron  gauge  is  a 
standard  thermocouple  gauge  with  one  probe  in  the  driven  tube 
and  another  monitoring  the  pressure  at  the  vacuum  pump.  Neither 
gauge  system  is  electrically  damaged  by  operating  at  atmospheric 
pressure.  So  far,  it  has  not  been  necessary  to  protect  any 
gauge  systems  by  isolation  during  firing.  This  is  important 
from  a  safety  standpoint-  as  well  as  an  aid  in  avoiding  possible 
readout  errors. 
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A- 5  VACUUM  SYSTEM 

A  satisfactory  range  of  initial  pressures  in  the  driven 
tube  was  deemed  to  be  from  50  microns  to  10  mm  Hg.  To  avoid 
the  additional  complication  of  a  diffusion  pump,  a  high  quality, 
high  capacity  mechanical  nump  is  used.  Though  pumping  speed  is 
slow  near  ultimate  pressure,  this  range  is  seldom  used  for 
decontaminating  the  tube  during  firing.  With  the  tube  clean 
enough  to  avoid  excessive  gassing  in  the  operating  pressure 
ranges,  the  simplest  and  quickest  way  to  achieve  a  contamina¬ 
tion-free  test  gas  is  by  repeated  filling  with  test  gas  and 
evacuating.  For  instance,  with  the  tube  at  atmospheric  and 
full  of  hydrogen  after  firing,  only  two  or  three  minutes’  pump¬ 
ing  will  evacuate  it  to  .1  mm  Hg.  By  refilling  with  test  gas 
and  evacuating  to  test  pressure,  remaining  hydrogen  contamina¬ 
tion  is  approximately  only  one  part  in  ten  thousand. 

With  the  addition  of  a  liquid  nitrogen  trap,  the  gassing 
contaminants  important  at  low  pressures  (including  pump  oil)  are 
removed.  However,  if  the  trap  is  used,  a  bypass  system  should 
also  be  employed  since  when  pumping  at  high  pressures,  severe 
and  wasteful  boiling  occurs  in  the  trap.  With  the  system  tight, 
pumping  for  twelve  hours  with  only  the  mechanical  pump  will 
reduce  the  pressure  to  a  few  microns.  Using  the  trap,  one 


micron  can  be  obtained  in  two  hours. 
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A- 6  INSTRUMENTATION 

The  instrumental ion  discussed  below  are  those  devices  and 
recording  instruments  pertinent  to  the  general  operation  of  the 
shock  tube.  The  microwave  instrumentation  is  discussed  in 
Section  A. 7.  The  specific  problem  under  study  requires  a  high 
degree  of  accuracy  in  all  measuring  and  recording  instruments. 
Three  specific  areas  will  be  discussed;  accurate  location  of 
the  aerodynamic  shock,  measurement  of  dynamic  pressures  and 
recording  of  fast  rise,  small  amplitude  electrical  signals. 

Location  of  Shock  Front 

Accurate  estimates  of  the  equilibrium  properties  of  the 
gas  behind  the  incident  shock  can  be  made  if  the  initial  pres¬ 
sure  and  the  shock  velocity  are  known,  the  latter  parameter 
being  the  most  sensitive.  Measurement  of  the  shock  velocity 
can  be  made  by  measuring  the  time  interval  between  two  points 
along  the  tube  where  an  observable  property  associated  with  the 
shock  front  can  be  detected.  The  difficulty  in  this  problem  is 
the  measurement  of  a  property  directly  attributed  to  the  shock 
front.  Basically  only  three  fundamental  parameters  are  im¬ 
plicitly  associated  with  the  shock;  a  temperature  change,  a 
pressure  change  or  a  density  change.  Measurement  of  parameters 
induced,  such  as  electron  density  and  light  radiation  to  locate 
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the  shock  front,  assumes  that  these  properties  abruptly  change 
across  the  front.  This  assumption  is  never  strictly  valid  but 
may  be  satisfactory  at  high  Mach  numbers,  when  the  use  of  ioniza¬ 
tion  gauges  or  photomultipliers  may  be  acceptable. 

If  it  is  desirable  to  detect  the  actual  aerodynamic  shock 
front,  one  of  the  fundamental  property  changes  must  be  observed. 

A  direct  measurement  of  the  density  change  can  be  made  optically 
by  measuring  the  deflection  of  a  light  beam.  However,  ever,  for 
moderate  strength  shocks,  this  method  is  very  delicate. 

A  direct  measurement  of  pressure  is  possible  except  that 
the  response  time  of  the  detector  strongly  limits  the  possible 
resolutions.  That  is,  using  a  good  quality  crystal  detector,  a 
rise  time  of  three  microseconds  is  typical.  If  the  shock  veloc¬ 
ity  is  on  the  order  of  3  mm/jisec  (~M9)  ,  the  spatial  resolution 
is  about  1  cm.  Assuming  the  rise  time  is  consistent  for  two 
gauges  and  the  shock  pressure  profile  is  constant  during  the 
interval  of  the  measurement,  a  'elocity  measurement  is  possible 
with  a  spatial  resolution  of  .5  cm.  Hence,  to  measure  velocity 
to  1%,  the  distance  between  gauges  should  be  .5  meters.  It  is 
clear  however  that  for  the  purpose  of  locating  the  shock  as  a 
method  of  correlating  events,  the  gauges  have  limitations.  In 
addition,  if  measurements  are  to  be  made  over  a  wide  range  of 
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conditions,  continual  adjustment  of  sensitivity  must  be  made. 

A  more  serious  difficulty  is  the  isolation  of  the  gauge  from 
the  mechanical  shock  generated  wher  the  diaphragm  breaks  (valve 
opens) .  If  the  mechanical  shock  reaches  the  gauge  prior  to  the 
arrival  to  the  aerodynamic  shock,  a  spurious  signal  may  render 
the  measurement  useless.  A  discussion  of  pressure  gauges  and 
measurements  will  be  made  in  a  following  subsection. 

Heat  Transfer  Gauges 

A  satisfactory  solution  to  a  detector  for  locating  the 
shock  front  has  been  found  to  be  a  heat  transfer  gauge.  The 
basic  gauge  is  a  thin  film  resistance  thermometer  element 
bonded  to  a  glass  slug.  The  gauge  response  is  a  function  of 
not  only  the  temperature,  but  also  the  density.  This  combined 
response  allows  the  useful  operation  from  M~4  to  ~  14  in  this 
facility.  The  gauge  is  fabricated  by  painting  a  thin  strip  of 
metallic  silver  coating  less  than  1  micron  thick  and  1  mm  wide 
on  a  glass  disk  \  inch  in  diameter.  This  film  is  then  baked, 
and  connecting  leads  are  fed  through  twTo  holes  in  the  glass  disc 
and  soldered  to  the  film.  The  'tire  front  surface  is  then 
given  a  thin  coating  of  magnesium  fluoride,  which  aids  in  pre¬ 
venting  the  collection  of  the  charge  diffused  in  front  of  the 
shock  (Ref.  1) .  Unsatisfactory  measurements  are  obtained  without 
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the  coating.  In  operation,  a  steady  current  of  10-50  ma  is 
passed  through  the  resistance  element  (thin  film) .  The  tempera¬ 
ture  change  as  the  shock  passes,  creates  a  change  in  resistance 
which  generates  a  change  in  voltage  if  the  current  is  nearly 
constant . 

An  extremely  simple  method  of  generating  an  essentially 
constant  current  is  to  use  a  large  voltage  with  a  large  series 
resistor  to  the  gauge  as  shown  in  Fig.  A-6.1.  To  be  sure,  this 
system  is  wasteful  of  power,  but  its  simplicity  over  a  suitable 
transistor  stabilized  constant  current  supply  for  each  gauge  is 
obvious. 


Fig.  A-6.1.  HTG  Power  Supply 


To  obtain  a  significant  pulse  (20  volts)  from  this  cir¬ 
cuit,  very  large  amplification  (20xl03)  of  the  voltage  pulse 
must  follow.  If  the  fast  response  of  the  gauge  is  to  be  fully 
utilized;  the  amplifier  should  have  a  rise  time  of  .3  pL sec . 
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To  make  the  gauge  useful  as  a  standard  and  reliable  shock 
front  detector  the  circuit  shown  in  Fig.  A-6.2  has  been  designed 
and  built.  The  circuit  is  intended  to  generate  a  sharp  pulse 
and  no  attempt  is  made  to  preserve  the  long  term  (10  jusec  or 
more)  variations  occurring  at  the  gauge.  That  is,  all  the  time 
constants  are  made  conveniently  small.  The  grounded  grid  first 
stage  utilizes  the  low  impedance  character  of  the  gauge  to  simul¬ 
taneously  provided  a  very’ stable  first  amplifier  and  the  required 
bias  current  for  the  gauge  element.  The  following  two  stages 
of  amplification  are  conventional  RC  pentode  amplifiers  utiliz¬ 
ing  an  extremely  high  transconductance  tube.  The  output  cathode 
follower  allows  extended  cabling  without  loss  of  response  time 
which  is  typically  .3  microsec.  The  combination  of  series  con¬ 
densers  and  the  normally  closed  relay  in  the  final  output  pro¬ 
vides  a  simple  method  of  avoiding  spurious  trigger  signals  which 
occur  from  the  various  switching  operations  performed  in  the 
firing  sequence.  Because  of  poor  grounding  practices  in  other 
nearby  facilities,  a  common  failing  in  many  laboratories,  minor 
power  line  noises  generated  by  switching  transients  are  ampli¬ 
fied  by  the  high  gain  amplifiers.  To  minimize  this  problem, 
the  first  step  in  the  firing  sequence  disconnects  all  unneces¬ 
sary  ac  components  from  the  system  (fans,  heaters,  pumps,  etc). 
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The  first  step  also  applies  a  dc  voltage  to  the  time  delay  cir¬ 
cuit  of  the  output  relay  which  opens  a  few  seconds  later.  This 
system  is  not  a  replacement  for  proper  grounding  procedures  but 
rather  an  insurance  to  avoid  60  cycle  pickup  and  line  transients. 

Care  must  be  taken  in  the  general  circuit  construction 
because  of  the  high  gain.  The  general  circuit  construction 
utilizes  v  ctor  turret  type  sockets.  Separate  filament  sup¬ 
plies  are  used  as  a  precaution  against  high  frequency  oscilla¬ 
tion.  The  separate  transformers  are  practical  if  multiple  units 
are  to  be  constructed  on  the  same  chassis.  For  one  amplifier 
unit  one  supply  is  probably  adequate.  The  B+  supply  can  be 
unregulated  with  minor  filtering.  At  300  v  >lts,  the  current 
consumption  is  about  70  ma  per  amplifier.  No  gain  control  is 
used  since  in  general  the  purpose  is  to  provide  a  fast,  large 
trigger.  For  most  of  the  shock  conditions  shewn  in  Fig.  B-4.5, 
the  amplifier  ir  saturated,  thus  generating  a  60  volt,  .2  jisec 
rise  time  pu?  se,  independent  of  shock  conditions. 

With  this  heat  transfer  gauge  and  amplifier,  spatial 
resolution  of  the  location  of  the  front  is  comparable  to  the 
actual  discontinuity  in  physical  parameters  occurring  across 
the  shock.  Typically,  this  resolution  is  1  mm,  indicating  that 
for  a  shock  velocity  of  3  mm/^jsec  the  required  distance  between 
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detectors  in  order  to  determine  velocity  to  1%  is  about  10  cm, 
.uming  the  timer  can  resolve  .3  ^sec  or  better. 

Pressure  Detectors 

Measurement  of  pressures  or  pressure  variations  behind  the 
shock  is  extremely  difficult  if  accuracy  and  space  resolution 
are  required.  The  only  suitable  gauge  elements  for  this  purpose 
are  piezoelectric  crystals.  The  base  problem  in  the  design  of 
a  gauge  for  shock  tuba  use  is  that  the  requirements  of  high 
natural  resonant  frequencies  and  high  sensitivities  are  con¬ 
tradictory.  The  second  problem  is  the  design  of  the  crystal 
mounting  in  the  gauge  so  that  the  output  is  insensitive  to  ex¬ 
ternal  str ac.ses  at  right  angles  to  the  measured  stress. 

If  one  considers  that  the  useful  pressure  information  in 
th:?  shock  tube  occurs  for  an  interval  of  a  few  hundred  micro¬ 
seconds,  there  is  a  partial  elimination  of  the  measuring  prob¬ 
lems  if  a  gauge  intended  only  for  dynamic  purposes  is  used. 

A  precision  calibrated  crystal  microphone  with  a  natural 
resonant  frequency  of  about  200  kc  is  the  Massa  model  M-213. 

This  microphone  has  good  sensitivity  and  directional  resolution, 
and  can  withstand  severe  overloads.  A  typical  pressure  trace 
using  this  microphone  is  shown  in  Fig.  A-6.3.  A  low  pass  filter 
of  RC  constant  20  ysec  has  been  used  to  remove  the  initial  ring¬ 
ing  of  the  crystal,  as  shown  in  Fig.  A-6.4, 
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Upper  Trace:  Massa  Microphone 
with  no  filtering  -  gain  1  5  db 
5  volt/ cm . 

Lower  Trace:  Electric  Field 
Detector  -  50  mv/cm. 

Sweep  Time:  20  (is/cm. 

Initir.j.  Pressure:  3  mm  Hg  -  Air. 
Shock  Velocity:  3000  meters/sec. 


Figure  A-6. 3 


Upper  T  race:  Heat  transfer  gauge 
current  10  ma,  5  mv/cm. 

Lower  Trace:  Massa  Microphone 
with  RC  =20  lis  filter  -  gain  20  db 
5  volts/cm. 

Sweep  Time:  50  jUs/cm. 

Initial  Pressure:  2  mm  Hg  - 


Figure  A-6. 4 


Shock  Tube  Instrumentation,  Typical  Traces 
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Since  this  crystal  microphone  is  only  us  ul  for  dynamic 
pressure  measurements,  its  calibration  is  difficult.  In  general 
the  calibration  supplied  with  the  microphone  is  adequate.  If 
doubtful,  the  microphone  can  be  mounted  in  an  inflated  balloon, 
the  pressure  statically  measured  with  respect  to  atmosphere,  and 
the  balloon  punctured.  The  resultant  step  pressure  change  sup¬ 
plies  a  positive  calibration  point.  Since  the  gauge  and  ampli¬ 
fier  are  relatively  inexpensive,  their  use  is  strongly  recom¬ 
mended  if  similar  dynamic  measurements  are  required. 
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A- 7  MICROWAVE  INSTRUMENTATION 

Microwave  Circuit 

The  circuit  is  as  shown  in  Fig.  A-7.1  includes  a  micro- 
wave  signal  generator  (2-4  kmc,  60  mw) ,  a  coupling  device 
delivering  this  signal  to  the  shock  tube,  detectors  of  electric 
and  magnetic  fields,  a  recording  system,  consisting  of  oscillo¬ 
scopes  and  photographic  cameras,  marker  generators,  wave  shapers 
and  delay  generators.  Sometimes  a  counter  is  used  to  measure 
the  time  between  markers  and  consequently  the  speed  of  the 
shock . 

The  microwave  signal  generator  is  a  reflex  klystron  with 
external  tunable  cavity.  This  klystron  is  powered  by  a  power 
supply  which  has  the  possibility  to  change  reflector,  grid  and 
beam  voltage.  The  frequency  stability  of  this  generator  was 
good  enough  for  our  purpose,  but  the  measurements  being  made 
were  particularly  sensitive  to  changes  in  the  amplitude  of  the 
signal.  Normally  the  klystron  is  switched  on  and  off  and  the 
maximum  tolerated  short  term  change  in  the  power  output  should 
be  below  2%.  Among  the  various  methods  tested,  the  one  using 
negative  pulses  in  the  grid  together  with  a  Zener-stabilized 
bias  power  supply  was  found  to  be  the  most  convenient. 


Fig.  A' 7.  1  Microwave  and  recording  c 
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Coupling  System 

The  output  of  the  klystron  is  fed  through  an  attenuator  to 
a  coaxial  cable  an  isolator,  and  from  this  to  the  coaxial  wave¬ 
guide  transition . The  shock  tube  is  equivalent  to  a  circular  wave¬ 
guide  operating  in  the  TEX 1  mode.  The  coaxial  to  waveguide 
transition  is  shown  in  Fig.  A-7.2.  Matching  is  made  by  chang¬ 
ing  the  penetration  of  the  antenna  and  the  position  of  the 
short  circuit  piston.  The  additional  short  circuit  stub  allows 
for  adjustment  of  the  bandwidth  changing  its  length.  The  shock- 
tube  must  be  vacuum  tight,  and  this  is  achieved  by  using  a 
phenolic  tube  around  the  antenna  and  a  bellows  on  the  piston. 

Matching 

The  matching  of  the  coaxial  line  to  waveguide  is  very 
important.  All  the  measurements  are  based  on  the  fact  that  the 
amplitude  of  the  microwave  signal  coupled  into  the  shock  tube 
is  constant.  This  condition  is  difficult  to  obtain  because  the 
plasma  produces  a  strong  reflection;  standing  wave  ratio  of 
twenty  is  not  unusual.  This  reflected  signal  must  be  absorbed 
without  reflection  in  the  isolator.  If  reflections  take  place, 
they  superimpose  on  the  incoming  signal  with  a  phase  that 
changes  as  the  plasma  moves.  This  effect  is  equivalent  to  a 
fluctuation  of  the  input  signal.  The  standing  wave  pattern 
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then  becomes  a  function  of  the  position  of  the  plasma  or  reflect-- 
ing  surface  which  is  undesired. 

Electromagnetic  detectors  mounted  on  the  waveguide  surface 
pick  up  the  sum  of  both  standing  wave  and  fluctuation  patterns 
both  having  the  same  periodicity.  A  detector  moving  with  the 
reflecting  surface  picks  up  the  fluctuation  only.  The  phase 
between  the  two  patterns  may  take  any  value:  if  it  is  coinci¬ 
dent  with  one.  for  instance  E ,  the  output  from  this  detector 
will  give  a  higher  standing  wave  ratio  chan  the  output  from  the 
other.  If  the  phase  angle  between  the  two  patterns  has  any 
value  the  resultant  pattern  is  not  symmetric.  In  the  computa¬ 
tion  of  the  plasma  properties  it  is  assumed  that  the  electro¬ 
magnetic  energy  is  perfectly  coupled  to  the  waveguide,  hence 
perfect  symmetry  must  be  obtained  to  utilize  the  measurement. 

Several  methods  have  been  used  to  get  the  desired  match¬ 
ing  but  the  only  one  accurate  and  simple  enough  is  the  follow¬ 
ing.  A  short  circuit  piston  is  displaced  along  the  shock  tube, 
simulating  a  strong  reflecting  plasma.  Moving  with  the  piston 
is  a  detector  picking  up  the  magnetic  field  near  the  wall  where 
it  is  maximum.  The  output  from  this  detector  gives  the  fluctua¬ 
tion  produced  by  the  mismatch  only.  Successive  adjustments  of 
the  antenna  to  waveguide  transition  are  made  attempting  to 
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eliminate  this  fluctuation,  thus  giving  a  constant  output  when 
the  piston  moves. 

In  order  to  reduce  the  frequency  sensitivity  of  the  an¬ 
tenna,  the  signal  from  the  generator  is  altered  slightly.  The 
antenna  short  circuit  stub  is  then  adjusted  until  a  minimum 
change  occurs  in  the  output  of  the  piston  detector  during  this 
frequency  change. 

The  limit  of  this  system  of  matching  comes  mainly  from 
incomplete  contact  between  piston  and  wall.  A  typical  value 
of  the  fluctuation  after  the  matching  has  been  performed  is  2%, 
but  this  corresponds  to  a  smaller  value  when  the  plasma  is 
reflecting  the  signal  because  of  the  significantly  better 
reflecting  properties  of  the  metallic  piston. 

Detectors 

The  electromagnetic  wave  propagated  inside  the  tube  has  a 
wavelength  X^  which  is  a  function  of  the  frequency,  mode  of 
propagation  and  diameter  of  the  tube.  Some  of  the  radiating 
energy  is  reflected  from  the  plasma  and  the  ratio  of  reflected 
to  incident  power  is  a  function  of  electron  density,  collision 
frequency  and  profile  of  the  plasma.  The  combination  of  inci¬ 
dent  and  reflected  waves  combines  to  form  a  system  of  standing 
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waves  along  the  tube.  This  pattern  has  a  period  of  . 


When 
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the  plasma  moves  the  pattern  also  moves  with  the  same  speed, 
provided  the  properties  of  the  plasma  remain  constant. 

A  device  detecting  the  electromagnetic  field  and  located 
in  a  fixed  position  detects  a  signal  whose  amplitude  goes  from 
a  maximum  value  equal  to  the  sum  of  incident  and  reflected 
fields  to  a  minimum  equal  to  the  difference  between  these  two. 

2v 

For  a  plasma  velocity  v  this  fluctuation  has  a  frequency  f  =  - —  . 

m  xg 

Two  types  of  detectors  have  been  used,  one  for  the  elec¬ 
tric  and  the  other  for  the  magnetic  field  amplitude,  as  shown 
in  Fig.  A-7.3.  The  electric  field  detector  is  located  in  the 
region  where  it  takes  its  maximum  value,  picking  up  the  com¬ 
ponent  parallel  to  the  antenna  and  using  a  small,  flush  mounted 
surface  in  order  to  avoid  introducing  anything  inside  the  shock 
tube  which  would  disturb  the  flow.  This  surface  is  connected  to 
a  coaxial  resonant  circuit  that  matches  impedances  to  the  crystal 
detector  and  short  circuits  the  currents  produced  by  the  flow  of 
charges  from  the  plasma  to  the  tube  walls.  The  magnetic  field 
detector  is  a  slot  excited  by  the  current  flowing  in  the  wall. 
Because  of  the  relation  between  magnetic  field  and  current  this 
slot  is  equivalent  to  a  loop  located  in  a  plane  perpendicular  to 
the  slot.  The  field  to  be  checked  is  the  transversal  component 
of  H.  The  slot  should  be  oriented  in  a  transverse  section  of 
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Figure  A-7.3  (a) 
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the  tube  and  at  the  po^nt  where  the  electric  field  is  maximum. 

The  signal  induced  in  the  slot  is  coupled  to  a  two  line  resonant 
circuit  c.nd  to  a  crystal  detector. 

For  both  E  and  H  crystal  detectors  an  output  is  obtained 
which  is  related  to  the  amplitude  of  the  respective  field.  The 
characteristics  of  the  crystals  cannot  be  considered  linear  nor 
quadratic  and  the  required  accuracy  demands  a  calibration.  The 
time  response  of  the  detectors  is  less  than  10“  sec  if  their 
output  is  directly  connected  to  cathode  followers.  For  simplic¬ 
ity's  sake  ac  amplifiers  are  used  exclusively  hence  the  dc  level 
is  not  present  However  switching  off  the  klystron  gives  this 
zero  level  reference.,  permitting  the  absolute  amplitude  deter¬ 
mination  from  the  single  trace.  The  markers  ara  equally  spaced 
at  known  intervals  providing  a  timing  reference  and  a  means  of 
compensating  for  ronlj nearities  in  the  sweep. 

To  perform  the  calculation  of  thv_  plasma  properties  the 
electric  and  magnetic  fields  should  be  determined  in  the  same 
transversal  section  of  the  tube.  A  special  section  has  been 
constructed  allowing  the  placing  of  the  two  detectors  plus  a 
shock  front  detector  in  the  same  transversal  section.  In  this 
arrangement  when  the  signal  from  one  detector  is  maximum  the 
other  is  minimum  and  vice  versa.  The  oscilloscopes  are  trig¬ 
gered  using  various  shock  front  detectors  as  discussed  in  A-6. 


41 


A- 8  DATA  RECORDING 

The  only  economically  feasible  method  of  recording  analog 
data  consistent  with  shock  tube  time  is  the  oscilloscope, 
although  it  : s  possible  to  record  some  data  with  video  type 
tape  recorders.  However,  the  cost  and  resolution  favor  the 
oscilloscope  and  camera.  While  this  basic  technique  is  common¬ 
place,  to  achieve  optimum  fidelity  of  the  data  consistent  with 
the  errors  in  the  system,  each  step  of  the  process  must  be  care¬ 
fully  analyzed.  This  error  analysis  is  in  many  instances  a 
function  of  the  recording  time.  The  range  of  interest  for  this 
facility  are  sweep  times  of  from  5  jisec  to  500  jisec. 

Two  types  of  distortion  can  be  delineated  in  the  over-all 
process;  these  are  amplitude  and  time  distortion.  To  trace  the 
source  of  these  and  estimate  their  relative  and  absolute  magni¬ 
tudes,  the  over-all  recording  process  is  divided  into  the  follow¬ 
ing  categories:  the  input  amplifier  up  t)  the  vertical  deflec¬ 
tion  plates  of  the  CRT,  the  sweep  generator  and  horizontal 
amplifier  up  to  the  horizontal  deflection  plates,  and  the  dis¬ 
tortion  in  the  CRT.  Also,  the  photographic  transfer  of  the 
CRT  face,  and  the  subsequent  extraction  of  numerical  data  from 
the  photographic  traces  must  be  considered.  By  far,  the  great¬ 
est  offender  in  these  distortion  processes  is  that  due  to  the 
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CRT  itself.  While  the  other  factors  are  not  always  unimportant, 
each  of  the  other  parts  of  the  process  except  numerical  data 
reduction  can  be  ma.i  stained  to  an  order  of  magnitude  better 
than  the  CRT  distortion  if  care  is  taken.  Even  with  the  best 
selected  tubes,  m  order  to  obtain  distortion  levels  on  the 
order  of  .5%.  only  the  central  region  of  the  CRT  screen  can  be 
utilized.  Of  course  by  so  doing,  this  also  tends  to  reduce  all 
the  other  distortion  causes  listed.  However,  if  these  processes 
are  independently  investigated  it  becomes  clear  that  the  CRT  is 
the  arch  criminal. 

The  specific  procedure  used  to  obtain  the  data  for  the 
shock  electron  density  profile  analysis  discussed  in  Section  B-3 
is  outlined  below.  The  critical  traces  were  made  using  selected 
CRT's  so  as  to  minimize  distortion.  The  photographic  transfer 
of  the  trace  was  made  using  a  modified  scope  camera  to  reduce 
the  image  size.  Polaroid  10-000  speed  film,  extremely  fast,- 
enables  the  use  of  a  small  lens  aperture.  With  this  arrange¬ 
ment  no  distortion  of  a  precision  grid  placed  on  the  CRT  was 
observable.  The  picture  obtained  was  concerted  to  a  transpar¬ 
ency  by  using  a  very  small  aperture  on  the  enlarging  lens.  The 
subsequent  transparency  was  enlarged  in  a  3"  x  4"  slide  projec¬ 
tor.  In  this  projector  a  3/8"  aperture  was  placed  inside  the 
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standard  lens  assembly  and  a  green  filter  was  added  to  reduce 
chromatic  abberation.  The  resultant  enlargement  of  about  12 
diameters  was  negligibly  distorted.  The  projection  screen  con¬ 
sisted  of  a  layer  of  fine  frosted  glass  backed  by  Plexiglass. 
The  Plexiglass  was  precision  scribed  with  .1"  squares.  The 
scribed  lines  were  painted  with  phosphorescent  paint.  By  using 
a  "black  lite"  fluorescent  behind  the  screen,  the  scribed  marks 
were  clearly  visible  and  direct  transfer  of  numerical  data 
from  the  picture  could  be  taken  by  working  from  behind  the 
screen.  Such  traces  have  also  been  made  on  large  graph  paper, 
but  such  paper  is  very  oftei.  distorted,  introducing  significant 
errors.  Under  the  best  conditions,  the  over-all  error  from 
amplifier  input  to  data  output  was  estimated  to  be  .4%  for  the 
maximum  values  increasing  to  2%  for  values  1/20  of  maximum. 

The  increase  in  error  becomes  significantly  large  for  small 
numerical  values  because  the  trace  cannot  be  resolved. 


CHAPTER  III 


ELECTRON  DENSITY  MEASUREMENT 


B-l  ANALYTICAL  FORMULATIONS  FOR  DETERMINING  ELECTRON  DENSITY 
The  following  is  a  summary  of  the  analytical  expressions 
used  to  determine  the  electron  density  in  the  combined  shock 
tube  -  waveguide .  The  polarized  radiation  is  introduced  in  the 
downstream  end  of  the  tube.  Probes  to  detect  the  radial  elec¬ 
tric  field  and  the  tangential  magnetic  field  are  upstream  of 
the  antenna  as  described  in  Section  A-7 .  To  determine  the 
pomt-by-point  variations  m  electron  density  as  the  plasma 
moves  past  the  detector  station,  expressions  are  required  relat¬ 
ing  the  measured  field  quantities  to  the  local  electron  density. 

For  the  TE- x  mode  of  propagation,  the  radial  electric 
field  and  the  tangential  magnetic  field  at  the  surface  of  the 
waveguide  are  given  by. 
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where  y  for  the  TE  ;  mode  is  1.842/R,  R  being  the  tube  radius. 
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The  function  u(x)  is  an  undefined  function  which,  however, 
contains  the  plasma  properties.  It  must  be  a  solution  to 
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where  primes  indicate  space  derivatives,  and  co  is  the  radiating 
frequency,  with 
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co  and  v  being  the  plasma  frequency  and  collision  frequency 
P 

respectively.  Assuming  the  field  components  are  measured  at 
right  angles,  with  tn<_  E^_  field  in  the  direction  of  polariza¬ 
tion  of  the  radiation,  then 

E  =  A  u  (x)  e1CJt  (B-5) 
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Since  the  values  of  the  fields  that  are  measured,  assuming  a 
square  law  detection,  are  the  square  of  the  real  parts  of  the 
respective  field  the  measured  quantities  are. 
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where  *  indicates  the  complex  conjugate.  Solving  for  K2  in 


Eq.  (B-3) 
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The  solution  of  K2  in  terms  of  the  field  values  of  Eqs 
'B*-'7)  and  (B-Cv  can  be  shown  to  be : 
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The  plasma  properties  m  terms  of  Re  K2  and  Im  K2  are  given  as. 
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The  values  of  the  field  in  the  above  equations  are  in  MKS 


units.  However  it  may  be  noted  that  each  field  quantity 
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appears  in  a  ratio.  All  that  is  required  for  solution  in  terms 
of  the  nondimensional  measured  quantities  are  the  appropriate 
scaling  ratios. 

If  the  measured  values  of  the  field  are  designated  as  E2 
and  H2 ,  then  the  distance  axis  scaling  needed  for  the  ratios 
-'3  7  and  -  ~^r~  can  be  determined  from  the  ratio  of  the  meas¬ 
ured  wavelength  Xg  (in  units  of  the  elongated  trace) .  That  is, 
if  the  horizontal  (distance  axis)  coordinate  is  measured  in 
units  of  Z  and  the  real  coordinate  in  X  units,  then. 
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where  X^  is  the  actual  waveguide  wavelength  (MKS  units) .  Hence 
to  scale  the  absolute  ratios  of  field  values  to  the  measured 
quantities. 
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To  scale  the  combined  ratios  of  electric  and  magnetic 
fields,  it  may  be  noted  that  in  the  free  space  region,  the 
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ratios  of  the  maximum  field  values  'occurring  180r  apart)  is 


given  by. 
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If  the  ratio  of  the  same  measured  quantities  is  designated  by 
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Then  Eqs.  (3-lCi  and  (8-111  can  be  put  in  terms  of  measured  quanti¬ 
ties  to  give 
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Hence  the  plasma  properties  can  be  obtained  through  the  use  of 


Eqs.  rr-12)  and  .‘B-13)  . 
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In  summary,  to  calculate  the  plasma  properties  at  any 
point  behind  the  shock  the  following  parameters  are  required. 


a  = 


H3 

max 

E3 

max 


Xg 

s  = 


x  = 


1.842 


a )  =  2nf 


waveguide  wavelength 

taken  in  the  free  space  region 

distance  scaling  ratios 

(X  measured  wavelength  on  enlarged  trace) 

g 

R  is  tube  radius 
f  radiating  frequency 


and  E2 ,  (E3)',  (E2)",  H2  ,  (H2 ) ' ,  the  field  values  and  derivatives 

at  the  wall  of  the  waveguide. 
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B-2  DATA  PROCESSING  PROGRAM 

In  the  preceding  section  describing  the  analytical  formula¬ 
tion  necessary  for  the  electron  density  measurement  the  basic 
quantities  required  for  solution  are  the  real  value  of  electric 
field  squared  (EE*)  its  first  and  second  derivatives,  t}ie  rial 
value  of  magnetic  field  squared  (HH*) ,  and  its  first  derivative. 

It  was  shown  that  no  absolute  measurement  of  these  quantities 

EE* 

was  required,  if  the  ratio  — - —  could  be  obtained  in  the  free 

HHmax 

space  region  in  front  of  the  shock.  The  precision  of  the  final 
calculation  is  extremely  sensitive  to  these  measured  quantities. 
While  it  has  been  indicated  in  Section  A-8  that  the  relative 
accuracy  for  part  of  the  recorded  data  is  about  .4%  it  is  im¬ 
portant  to  note  that  this  accuracy  is  related  to  the  data  taken 
at  that  instant  of  time.  Because  all  of  the  components  of  the 
measuring  system  are  constantly  varying,  it  is  essential  that 
all  of  the  critical  information  be  obtained  from  one  set  of 
traces  taken  at  one  time. 

The  greatest  variables  m  this  system  are  the  crystal 
detectors,  although  sweep  rates  and  amplifier  gains  may  also 
change.  If  the  desired  signal  to  be  detected  is  for  instance 
EE*,  the  crystal  detector  must  have  a  square  law  behavior. 

Since  this  requirement  is  never  precisely  true  a  correction 
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must  be  made  on  the  output  signal.  It  will  be  shown  how  this 
correction  can  be  made  directly  from  the  numerical  data  derived 
from  the  single  trace,  thus  avoiding  the  need  for  corroborating 
data  taken  at  a  different  time. 

The  calculation  of  electron  density  is  divided  into  three 
operations.  First,  the  raw  numerical  data  taken  from  the  scope 
trace  is  corrected  for  the  non-square  law  behavior  of  the 
crystal.  Second,  from  this  corrected  data,  first  and  second 
derivatives  are  computed.  Third  the  corrected  data  and  deriva¬ 
tives  are  combined  point“by-point  as  prescribed  by  the  analytic 
equations  of  the  preceding  section. 

Processing  of  the  Raw  Data 

The  basic  input  for  the  calculation  is  approximately  150 
data  points  from  each  of  the  raw  electric  and  magnetic  field 
detector  output  traces.  These  traces  include  one  full  period 
of  the  sinusoidal-like  function  preceding  the  plasma. 

Because  careless  errors  may  be  introduced  in  the  trar scrip 
tion  or  punching  of  the  raw  data,  the  program  first  checks  if 
successive  points  arc  monotonically  increasing  or  decreasing. 

The  object  of  this  procedure  is  to  identify  data  that  is  clearly 
impossible,  and  allow  corrections  or  changes. 
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Once  the  data  passes  this  initial  test  the  points  are 
smoothed  by  a  least  square  fitting  procedure  In  this  computa¬ 
tion  eight  points  are  used  to  obtain  the  constants  of  a  second 
power  equation  by  least  squares  fit.  Having  computed  the  equa¬ 
tion  constants  the  values  of  the  function  for  f  center  two 
points  are  calculated  and  stored  as  smooth  data.  The  program 
then  jumps  two  points  and  recomputes  the  least  square  fit  for 
the  next  eight  points  which  of  course  includes  six  from  the  pre¬ 
vious  calculation.  After  this  procedure  has  been  performed  for 
all  data  points  the  initial  raw  data  is  destroyed  and  the 
smoothed  values  inserted  in  its  respective  place.  The  need 
fcr  smoothing  the  data  is  comparable  to  the  use  of  a  "French" 
curve  m  graph  plotting.  The  basic  assumption  is  that  the  data 
is  originally  smooth  and  that  conversion  to  digits  has  intro¬ 
duced  a  random  error  in  the  least  significant  place.  By  reduc¬ 
ing  these  errors  an  improvement  can  be  obtained  m  the  subse¬ 
quent  derivatives  that  must  be  taken. 

The  next  step  is  to  modify  the  data  for  the  incorrect 
response  of  the  detector.  This  correction  procedure  is  based 
on  the  fact  that  the  data  function  outside  the  plasma  must  be 
a  sinusoid  i  '  the  detectors  are  measuring  EE1,  The  basic  prob¬ 
lem  is  to  find  a  multipli ca  :: *-e  function  [g,x)j  such  that  when 
multiplied  by  the  data  furcnioo  i  f 'x. j  the  resulting  function 
will  be  a  sinusoid  ele^  if ed  ebo^e  zero 
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i.e.,  g(x)  f(x)  =  A  sin  0  +  A  +  a 

The  difficulty  in  determining  the  g(x)  function  is  that  since  a, 
which  is  related  to  the  standing  wave  ratio,  is  unknown,,  then 
there  are  an  infinite  set  of  g(x)  functions  which  will  satisfy 
the  above  criteria.  The  solution  to  this  dilemna  is  found  in 
the  fact  that  the  crystal  nonlinearity  (to  the  square  law)  must 
be  a  smooth  and  slowly  varying  function,  a  fact  clear  from  the 
physics  of  tb,_  semiconductor.  If  a  is  on  the  order  of  A/3  or 
smaller,  the  latter  condition  becomes  a  very  stringent  limit  on 
possible  g(x)  functions.  The  first  step  to  determine  the  g(x) 
function  is  to  determine  the  period  of  the  sinusoidal-like  func¬ 
tion  f(x)  in  the  free  space  region.  This  is  accomplished  by 
computing  a  function  like  a  derivative  over  the  entire  curve. 
This  derivative  procedure  consists  of  computing  multiple  dif¬ 
ferences  at  each  point  and  averaging  the  result.  The  reason 
for  this  technique  is  to  establish  the  zero  slope  points  with 
the  least  error  being  introduced  from  the  random  data  scatter. 

After  locating  the  first  two  zeros,  forty  points  are  gen¬ 
erated  by  interpolation  between  zeros,  which  of  course  corre¬ 
spond  to  the  first  maximum  and  minimum  or  the  reverse  of  the 
original  curve.  The  forty  points  are  chosen  such  that  if  the 


function  were  a  true  sine,  then  the  increment-  of  the  function 
between  points  would  he  equal.  This  is  done  so  that  the  g (x) 
function  will  not  he  unduly  weighted  by  a  preponderance  of 
points  at  each  end.  The  first  trial  g(x)  function  is  generated 
by  assuming  the  perfect  sine  curve  has  the  same  magnitude  as 
the  data  function.  Hence  the  gfx)  function  at  maximum  and  mini 
mum  ( a  and  B  in  Fig.  B-2.1)  is  identically  one.  This  then  pro¬ 
duces  a  g(x)  typically  like  that  in  Fig.  B-2.2.  On  the  g(x) 
curve,  the  region  from  A  to  D,  where  D  is  arbitrarily  1/2  f(x), 
is  least  square -fitted  by  a  second  order  polynomial.  Point  C, 
the  value  of  g(x)  for  f(x)  min  is  computed  from  the  curve  fit 
constants  and  this  value  is  used  to  compute  the  value  of  a. 


Distorted  Sine  Function  g(x)  Determination 


F ig .  B-2 . 2 


55 


Using  this  value,  a  new  magnitude  A  is  computed  and  the 
entire  g(x)  function  for  each  point  of  f(x)  is  generated.  This 
second  trial  g{x)  function  is  again  curve -fit  to  a  second  order 
polynomial  by  method  of  least  squares.  Prom  the  constants  of 
the  polynomial,  the  value  of  g(o)  and  the  standing  wave  ratio 
are  calculated  and  the  difference  between  the  curve  fit  g(x)  and 
the  trial  computed  g(x)  points  is  printed.  At  this  time  in  the 
proceeding,  the  computer  is  tired  and  returns  control  to  the 
operator.  The  operator,  making  a  human  type  judgment  on  the 
basis  of  the  differences  in  the  last  points.,  may  choose  a  dif¬ 
ferent  value  by  a  few  percent  of  g(o)  and  the  above  process  will 
be  repeated. 

It  is  desirable  that  the  companion  curve,  HH*,  to  the 
EE*  curve  be  processed  up  to  the  determination  of  the  g(x)  cor¬ 
rection  function  as  indicated  above.  If  the  g(x)  functions  have 
been  properly  chosen  for  each  curve,  the  standing  wave  ratio 
will  be  identical.  With  a  little  practice  the  standing  wave 
ratios  can  be  made  to  agree  to  better  than  2%. 

Once  the  g(x)  functions  are  satisfactorily  determined,  all 
the  data  is  corrected  thus  yielding  the  "true"  values  of  EE*  and 
HH*  along  the  curves.  By  a  method  of  multiple  differences,  the 
derivative  of  these  functions  are  taken  yielding  (EE*)',  (EE*)" 


and  (HH*)  '  . 
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The  computer  program  thus  far  described  is  intended  to 
take  data  derived  from  a  single  trace,  apply  corrections  and 
produce  derivatives  at  fixed  intervals.  The  functions  gener¬ 
ated  must  not  only  be  self-consistent,  but  also  should  be  con¬ 
sistent  with  the  data  derived  from  the  companion  curve.  (EE* 
and  HH*  form  a  set.)  It  has  been  shown  how  the  program  "forces" 
an  over-all  amplitude  correction  by  assuming  the  deleted  func¬ 
tion  in  the  free  space  region  is  sinusoidal.  This  type  of  cor¬ 
rection  will  compensate  for  most  amplitude  distortions  due  to 
any  cause  that  is  not  a  function  of  the  time  coordinate  (hori¬ 
zontal  position  on  the  CRT) . 

Additional  information  which  can  be  obtained  from  the  data 
is  the  relative  wavelength  of  the  sinusoidal  function  in  the 
free  space  region.  It  may  be  recalled  that  the  wavelength 
determination  was  necessary  to  establish  the  period  of  the 
"forcing"  sinus  function.  Since  the  waveguide  wavelength  of 
the  microwave  radiation  can  be  accurately  measured  prior  to 
testing^  it  is  used  as  the  only  additional  piece  of  information 
that  is  not  directly  obtained  from  the  picture.  If  the  meas¬ 
ured  wavelength  and  the  relative  wavelength  are  combined,  an 
over-all  scaling  is  achieved.  In  addition,  if  the  scaling  fac¬ 
tors  for  the  EE*  and  HH*  are  compared,  an  accurate  contraction 
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or  expansion  factor  can  obtained  and  used  to  provide  exact 
time  correlation  between  points  on  either  curve.  If  one  addi¬ 
tional  marker  (such  as  the  location  of  the  shock  front)  is  com¬ 
mon  to  both  curves  an  exact  relation  between  the  curves  can  be 
made.  If  the  trace  also  contains  time  markers,  the  shock  veloc¬ 
ity  can  be  obtained.  It  may  be  noted  that  the  time  markers  also 
provide  the  location  of  the  zero  amplitude  axis  at  the  instant 
the  trace  was  made,  a  procedure  which  removes  the  uncertainty 
due  to  drift  in  all  systems. 

An  additional  piece  of  information  required  for  the  cal¬ 
culation  of  electron  density  is  the  ratio  of  the  maximum  on  the 

electric  field  trace  in  the  free  space  region  to  the  maximum  on 

EE*ax 

the  magnetic  detector  trace  -nnj —  .  This  ratio  provides  all  the 

“inax 

scaling  information  needed  to  use  the  nondimensional  curve  data 
in  an  absolute  determination  of  electron  density.  It  may  be 
recalled  that  the  maximum  values  required  were  obtained  to  pro¬ 
vide  amplitude  information  for  the  "forcing"  sinusoidal  correc¬ 
tion  function. 


Hence,  it  has  been  demonstrated  that  the  complete  set  of 
information  required  to  cc  pute  the  electron  density,  with  the 
exception  of  the  waveguide  wavelength  is  directly  obtainable 
from  the  two  detector  curves.  The  advantage  of  this  procedure 
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is  that  no  absolute  measurements  are  required <  with  the  excep¬ 
tion  of  the  waveguide  wavelength  and  interval  between  time 
markers.  (The  latter  is  used  only  for  the  velocity  determina¬ 
tion.)  This  clearly  implies  that  no  absolute  calibration  and 
maintenance  of  that  calibration  is  needed  for  the  instruments 
used  in  recording.  Except  for  the  requirement  of  a  linear  time 
scale  (assuming  shock  velocity  constant) ,  an  over-all  correction 
of  most  of  the  distortions  occurring  in  ail  phases  of  the  data 
collection  process  is  obtainable  from  the  data. 

All  of  the  informat  '■'n  obtain?  ■■  in  "data  proct.  •  Li  " 

phase  of  the  program  is  subsequently  used  in  the  point-by-point 
determination  of  electron  density  as  per  the  formu1 ations  of 
the  preceding  section.  Since  this  computation  is  entirely 
straightforward  no  discussion  will  be  given  here. 
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B-3  EXPERIMENTAL  DETERMINATION  OF  THE  ELECTRON  DENSITY 

BEHIND  A  NORMAL  AIR  SHOCK 

The  following  is  a  brief  description  of  the  experimental 
conditions  for  which  the  longitudinal  electron  density  profile 
has  been  determined.  A  series  of  six  independent  sets  of  data 
have  been  processed  corresponding  to  nearly  identical  chock  con¬ 
ditions.  The  purpose  of  the  nearly  redundant  data  was  to  check 
the  over-all  repeatability  and  the  effect  of  small  variations 
in  shock  condition  about  the  mean. 


The  specific  shock  conditions  are  tabulated  below. 


Case 

Shock  Velocity 
(meter/sec) 

Initial  Pressure 
(mm  Hg) 

Note 

1 

3100 

2.5 

Upper  curve 

2 

2900 

1.5 

Lower  curve 

3 

3030 

2.0 

4 

3015 

2.05 

Mean  values 

5 

3035 

1.95 

plotted 

6 

3020 

2.0 

The  velocities  quoted  above  were  calculated  directly  from  the 
electric  and  magnetic  field  detector  traces  as  previously  de¬ 
scribed  in  Section  B-2.  The  microwave  instrumentation  used  is 
describee  in  Section  A-7 .  The  specific  waveguide  w wo length 


for  all  of  the  runs  was  20.2  cm. 
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One  typical  set  of  traces  is  shown  in  Fig.  B-3.1.  Note 
the  periodic  time  markers  which  also  give  the  zero  level  of  the 
signals.  The  location  of  the  shock  front  was  determined  from  a 
third  trace. 

The  numerical  data  from  these  curves  was  obtained  in  the 
manner  described  in  Section  A-8.  Subsequent  numerical  analysis 
was  carried  out  using  the  computer  program  described  in  Sec¬ 
tion  B-2.  The  results  for  the  six  cases  are  summarized  in 
Fig.  B-3.2.  Here  the  mean  curve  represents  the  results  for 
cases  3  through  6  which  correspond  to  the  same  velocity  condi¬ 
tions  ± . 3% .  Cases  1  and  2  correspond  to  higher  and  lower  veloc¬ 
ities  respectively.  It  may  be  noted  that  the  behavior  of  the 
electron  density  is  in  correspondence  with  these  conditions. 

The  values  of  equilibrium  electron  density  were  calculated  for 
these  shock  conditions  using  the  procedure  outlined  in  B-4  and 
Fig.  B-4. 2. 

Because  the  shock  front  marker  detector  used  in  these 
runs  was  derived  from  a  pressure  detector  as  described  in  Sec¬ 
tion  A-6,  a  small  displacement  uncertainty  existed  in  the  re¬ 
sults  consistent  with  the  "jitter"  of  the  relatively  slow  rise 
of  the  detector.  The  "mean"  curve  in  Fig.  B-3.2  represents  a 
repositioning  of  the  separate  results,  where  the  maximum  shift 


corresponded  to  .75  cm. 
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Set  of  Electric  and  Magnetic  Field  traces  used  as  starting 
data  in  the  calculation  of  tho  electron  density  profile  behind  the 
shock.  (Case  4  in  table). 


Figure  B -3,  1 


THIS  FACE  INTENTIONALLY  LEFT  BLANK 
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he  test  gas  for  all  of  these  cases  was  commercial  com¬ 
pressed  breathing  air.  The  basic  result  has  not  been  signifi¬ 
cantly  modified  by  any  of  the  following  procedures. 

(1)  Addition  of  a  liquid  nitrogen  "freezeout"  trap  on  the 
test  gas.  Such  a  trap  would  of  course  remove  all 
water  and  o^l  present. 

(?)  Liquid  nitrogen  trap  in  the  vacuum  system  to  remove 
any  trace  of  pump  oil  from  the  tube. 

(3)  Mixing  of  high  purity  oxygen  and  nitrogen  in  propor¬ 
tion  to  "normal"  air  and  slight  variation  about  the 
standard  mixture.  These  gases  were  also  liquid 
nitrogen  trapped  to  remove  water  or  oil. 

(4)  The  design  of  the  mechanical  valve  discussed  in  Sec¬ 
tion  A-3  was  motivated  by  the  desire  to  remove  the 
possibility  of  cor  .  mination  due  to  the  use  of  plas¬ 
tic  diaphragms.  However,  no  significant  change  of 
the  electron  density  profile  was  observed  after  the 
conversion  to  the  mechanical  valve. 
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B-4  SUPPORTING  EXPERIMENTAL  AND  ANAI.YTICAL  INFORMATION 

Radial  Uniformity  of  the.  Plasma 

One  assumption  made  in  the  analysis  of  Section  B-l  for 
the  electron  density  was  that  of  radial  uniformity  of  the  plasma 
density.  Such  uniformity  can  be  altered  by  two  factors.  First 
is  the  presence  of  a  boundary  layer  along  the  walls  in  the  re¬ 
gion  behind  the  shock.  However,  for  the  penetration  distances 
indicated  by  the  results,  this  should  have  an  insignificant 
effect  on  the  assumption  of  uniformity. 

The  second  effect  is  the u.  of  ambipolar  diffusion.  A  cal¬ 
culation  of  this  effect  for  two  operating  conditions  has  been 
performed  as  summarized  in  the  following  subsection.  The  re¬ 
sult  indicates  the  assumption  of  radial  uniformity  is  basically 
correct . 

The  ambipolar  diffusion  equation  for  electrons  in  cylin¬ 
drical  coordinates  is: 


,  a  Bn 
Id  e 

—  ~r  r  — — 
r  dr  ,  or 


u  Bn  n 

n  _ e  _e 

D  Bz  D 


0 


(B-15) 


\llc  T 

where  D  =  — -  is  the  ambipolar  coefficient .No  4>  dependence  of 

Mi?'i 

the  electron  distribution  is  included. 


Only  diffusion  normal  to  the  wall  will  be  studied  and  for 


that  particul.  '  case  Eq.  (B-15)  becomes: 


65 


1  _d_ 
r  dr 


dn  \  n 
e  e 


0 


(B-16) 


The  term,  n  which  is  related  to  the  creation  of  electrons  due 
e 

to  chemical  reactions  inside  the  ionized  gas,  as  yet  must  be 
defined. 

The  mechanism  for  the  creation  of  the  electron  is  given  by 
the  motion  of  a  shock  wave  along  a  tube.  From  the  several  reac¬ 
tions  which  can  generate  electrons  in  air  due  to  t  effect  of 
the  shock  wave-  the  most  important  is: 


Therefore - 


N  +  0  -  NO  +  e 


n  =  k  (K  n„n  -n  2 ) 
e  r  c  N  0  e 


Replacing  (B-18)  into  (B-16)  we  have: 


(B-17) 


(B-18) 


<i  a  dn 
_1  _o_  _ e 

r  dr  \  dr 


(K  n  n  -n 
c  N  0 


0 


(B-19) 


Eq.  (B-19)  can  be  written  in  the  form: 


1  _d_ 

r 

dr 

k  n 
r  eo 

K=Vo  J 

r  dr 

+  D 

neQ2  I 

(B-20) 


eo 


where  n 


is  the  equilibrium  electron  number  density. 
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It  is  known  from  chemistry  considerations  (see  Penner, 
Chemistry  Problems  in  Jet  Propulsion) ,  that, 


K  n„n. 
c  N  0 


=  1 


n 


eo 


(B-21) 


Then  Eq.  (B-20)  becomes: 


1  J_ 

r  Sr 


r  & 
Sr 


k  n 
r  eo 


(1-y2 )  =  0 


(B-22) 


This  equation  has  a  boundary  condition  given  by  the  following 


relation  (see  Section  D-3) . 


D 


n 

e 


dn 

_ e 

dr 


(B-23) 


which  can  be  written, 


D  dy  =  /kT  (B_24) 

y  dr  */  m^ 


The  Eq.  (B-22)  with  the  condition  (B-24)  must  be  solved  numerically 
and  it  has  been  done  for  5  shock  wave  configurations.  The  data. 


corresponding  to  these  configurations  appear  in  Table  I. 


TABLE  I 


Shock 

Velocity 

(m/sec) 

- 

T 

(°K) 

P/Po 

neo 
(cm"3  ) 

Of 

(cm) 

2750 

2790 

8.45  x  10"3 

4.4  x  109 

1.404 

3000 

3000 

1.00  x  10"2 

2.10  x  101C 

.639 

3250 

3210 

1.08  x  10"3 

7.2  x  1010 

.348 

3500 

3390 

1.22  x  10-2 

2  x  1011 

.204 

3750 

_ _ 

3660 

1.30  x  10“2 

6.2  x  1011 

.120 

Of  primary  interest  is  the  effect  of  the  coefficient, 


a2 


k  n 
r  eo 


(B-25) 


upon  the  electron  distribution. 

The  coefficient  a  has  the  dimension  of  a  length  and  can 
be  considered  to  be  a  characteristic  diffusion  length.  Table  I 
shows  that  a  changes  greatly;  it  implying  that  for  the  cases 
where  a  is  small,  the  diffusion  will  be  very  slow  at  the  begin¬ 
ning  but  will  fall  off  sharply  near  the  wall  of  the  cylinder. 

This  does  not  happen  for  much  higher  values  of  a  (see  Pig.  B-4.1). 
The  experimental  investigation  made  at  GASL  for  a  case  where  the 
shock  wave  velocity  was  3100  m/sec  corresponds  to  large  values 
of  a.  For  this  velocity  case,  the  transverse  gradients  are 
small  over  most  of  the  tube  radius,  becoming  large  only  in  the 


thin  outer  sheath. 


Crn 
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Microwave  Interferometer 

A  conventional  phase  detection  interferometer  was  adapted 
to  a  special  test  section  of  the  shock  tube.  This  test  section 
was  made  of  thin  glass  with  the  identical  inside  diameter  of  the 
shock  tube.  A  supporting  housing  provided  the  necessary  struc¬ 
tural  capacity  and  housed  the  microwave  horns.  The  nominal  wave¬ 
length  of  the  radiation  was  about  1  cm.  The  nominal  beamwidth 
resolution  was  on  the  order  of  2.5  cm. 

With  this  instrument,  the  shock  tube  was  operated  using 
air.  For  the  range  of  pressures  and  velocities  indicated  in 
Fig.  B-4.2,  the  electron  density  was  within  20%  of  the  equilib¬ 
rium  values.  Further,  the  change  across  the  shock  front  with 
the  probable  2.5  cm  resolution  appeared  essentially  as  a  step. 

Sensitive  Check  on  the  Equilibrium  Electron  Density 

By  assuming  a  simple  model  for  the  plasma  behind  the  shock 
to  be  a  sharp  transition  from  free  space  to  equilibrium  electron 
density  levels,  a  formulation  for  the  standing  wave  ratio  due  to 
the  addition  of  the  transmitted  and  reflected  signals  from  such 
a  plasma  can  be  easily  made. 

By  combining  the  formulation  of  the  SWR  with  the  properties 
of  the  gas  behind  an  air  shock  as  a  function  of  shock  velocity 
and  initial  pressure,  a  comparison  with  experimental  results  can 
be  made.  This  procedure  is  summarized  below. 
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First,  the  standinq  wave  ratio  will  be  calculated  in  terms 
of  the  electron  density  and  collision  frequency  which  are  assumed 
uniform  in  the  plasma  region. 

The  generalized  propagation  constant  as  defined  by 


d2  E 
dx: 


a  +  K2E  =  0 


is  given  by, 


K_  I  ^ 

=  / 


of. 

a 


1- 


03 

3 

‘ 

0 

oP 

[l-i  ^ 

03 

-  V‘ 


where 


03  =  the  propagating  frequency 

ajp  =  the  plasma  frequency 
v  -  the  collision  frequency 

y  =  the  root  of  the  Bessel  function  for  the  TEj  l  circu¬ 
lator  waveguide  mode  divided  by  the  tube  radius  . 


The  propagation  constant  can  be  formed  as  follows.  Define 


where 


K  "  k  K 
o  g 


71 


Now  let  the  radial  electric  field  and  the  tangential  magnetic 
field  in  the  free  space  region  be  given  by 

-ik  x  -ik  x 

o  o 

E  =  A  e  -r  B  e 

r 


-k  x  ik  x 

o  _  o 
A  e  -Be 


and  in  the  uniform  plasma  region, 

ikxx 

E  =  C  e 
r 


Normalizing  A  =  1 ,  the  boundary  conditions  are  matched  yielding. 


with 


E 

r 


-ik  x 


e 


o 


+  0 


ik  x 
o 
e 


0 


1  -  K 

_ a 

1  +  K 

g 


Thus,  the  real  component  of  the  field  produced  is, 


-ik  x  ik  x 

E  E  *  =  1  +  00*  +  00*  e  °  e  ° 
r  r 
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This  is  reduced  to, 


EE* 
r  r 


1 


±  V  (Re0)2  +  (lm9)a 


The  two  solutions  give  the  maximum  and  minimum  values.  Thus  the 
standing  wave  ratio  is, 


SWR 


1  +  V  (R  6)2  +  (I  0)2 
_ e _  m 

1  =  V  (Re0)2  +  ( Im0 ) 2 


where  0  is  a  function  of  ui  and  v. 

P 

By  solving  the  normal  set  of  aerodynamic  equations  and  using 
a  programmed  Mollier  diagram,  the  conditions  behind  the  shock  can 
be  given  as  a  function  of  shock  velocity  and  initial  pressure. 

To  simplify  the  calculation  of  equilibrium  electron  density 
from  temperature  and  density,  standard  tables  were  programmed  for 
the  computer.  Thus  as  a  function  of  shock  velocity,  the  equilib¬ 
rium  electron  density  is  computed  with  the  initial  pressure  as  a 
parameter.  Figure  B-4 . 2  shows  this  plot.  The  collision  fre¬ 
quency  is  computed  from  the  approximation 


V 


5  x  1Q1SP 


with  p  the  pressure  in  atmospheres,  T  in  °K. 


74 


If  now  these  plasma  properties  are  used  to  compute  K^,  the 
standing  wave  ratio  can  be  calculated  using  the  previous  formula¬ 
tion  as  a  function  of  shock  velocity  with  the  initial  pressure 
as  a  parameter  (Fig.  B-4.3).  Hence  as  a  function  of  accurately 
determinable  quantities,  the  electromagnetic  characteristics  of 
the  entire  system  can  be  computed  and  subsequently  compared  with 
the  experimental  quantities. 

It  should  be  appreciated  that  such  a  comparison  is  a  very 
sensitive  check  on  the  assumptions  made.  First,  the  aerodynamic 
properties  behind  the  shock  must  precisely  follow  the  gas  dynamic 
theory.  Second,  the  electron  density  and  collision  frequency 
must,  on  the  average,  be  identical  with  the  equilibrium  values. 
Third,  the  assumed  geometric  uniformity  of  the  plasma  region 
must  be  nearly  correct.  Fourth,  the  calculation  of  the  SWR 
based  on  the  simple  step  model  should  be  right.  When  comparing 
with  experimental s,  the  variables  measured  must  be  identically 
the  sirae  as  those  used  in  the  formulation.  For  instance,  the 
electric  field  detector  after  correction  should  indeed  be  meas¬ 
uring  only  EE*. 

In  this  light,  the  correspondence  between  the  analytically 
determined  and  the  measured  SWR  as  presented  in  Fig.  B-4.4  is 


remarkable.  Note  here  that  since  the  driver  pressure  is  constant 
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for  all  these  tests,  in  order  to  change  the  shock  velocity  there 
is  a  corresponding  change  in  initial  pressure.  Hence,  Fig.  B-4.4 
corresponds  to  a  cross  plot  of  Fig.  B-4.3  using  the  operating 
characteristic  curve  in  Fig.  B-4.5.  This  agreement  provides 
probable  verification  of  the  simple  analytical  model  used  in  the 
computation.  The  basic  assumptions  of  that  model  are  uniform 
plasma  densities  both  radially  and  transversely  and  that  the 
density  levels  are  as  coir  ated  from  the  equilibrium  thermodynamic 


state  of  the  shocked  gas. 
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CHAPTER  IV 

COMPARISON  BETWEEN  EXPERIMENTAL  AND  ANALYTICALLY 
DETERMINED  ELECTRON  DENSITY  PROFILE 

C-l  INTRODUCTION 

Historically  the  shock  tube  has  been  one  of  the  primary 
tools  used  to  investigate  the  non-equilibrium  air  chemistry. 
Like  most  scientific  endeavors,  the  purpose  of  such  research  i 
to  provide  a  model  consistent  with  all  the  data  which  may  be 
used  to  predict  the  behavior  of  che  system  under  a  variety  of 
conditions.  The  obvious  complexity  of  the  air  system  and  the 
difficult  experimental  techniques  used  for  investigation  have 
severely  limited  the  extent  and  reliability  of  the  input  data. 

One  of  the  important  tests  of  the  air  chemistry  model  is 
the  prediction  of  the  non-equilibrium  electron  density.  On  the 
basis  of  this  test  disagreement  has  been  found  between  the 
experiments  conducted  herein  and  the  "state  of  the  art  model" 
that  has  beer,  devised. 

The  following  sections  detail  the  digital  computer  pro¬ 
gram  used  to  operate  on  the  air  chemistry  model.  Secondly  the 
air  chemistry  model  is  considered  with  the  addition  of  impuri¬ 
ties  to  ascertain  their  effects.  Thirdly,  the  experimental 
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spectroscopic  investigation  of  the  shock  produced  plasma  is  dis¬ 
cussed.  And  in  section  four,  the  results  of  the  experimental 
and  analytical  programs  are  compared.  Because  of  the  complexity 
of  the  problem  and  the  inclusive  data  obtained  from  all  sources, 
it  would  be  reckless  to  draw  firm  conclusions  from  the  material 
presented  here.  However  it  is  probable  that  the  observed  dis¬ 
crepancy  may  be  fully  acounted  for  -  in  terms  of  impurities. 
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C-2  ONE-DIMENSIONAL  NON-EQUILIBRIUM  AIR  CHEMISTRY  PROGRAM 

(REF,  2) 

As  an  integral  part  of  the  work  performed  in  this  investi¬ 
gation,  a  digital  computer  program  has  been  developed  to  study 
one-dimensional  chemically  reacting  flows.  The  characteristic 
features  of  this  program,  as  well  as  the  philosophy  behind  it, 
will  now  be  described. 

To  one  familiar  with  the  problem,  it  is  evident  that  the 
program  should  be  so  structured  that  the  number  and  kinds  of 
species  and  chemical  reactions  can  be  easily  changed.  The  chem¬ 
ical  rate  constants  should  also  be  amenable  to  a  quick  altera¬ 
tion  if  desired.  The  reason  for  this  stems  from  the  complexity 
of  the  chemistry  and  also  the  uncertainty  of  the  related  data. 
Much  of  the  empirical  data  is  obtained  from  extrapolations  of 
data  into  ranges  of  temperature  and  into  situations  which  differ 
appreciably  from  the  conditions  under  which  the  original  measure¬ 
ments  were  made.  Such  limitations  suggest  that  provision  be 
made  for  varying  or  correcting  the  data  in  a  routine  manner. 

Accordingly,  the  present  study  has  as  its  purpose  the  com¬ 
plete  formulation  of  the  problem  of  generalized  one-dimensional 
chemically  reacting  flows  with  molecular  vibrational  relaxa¬ 
tion,  and  includes  the  lowing  features: 
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Optional  Relaxation  of  Vibrational  Energy  of  Diatomic 
Species 

The  diatomic  molecules  of  a  given  species  are  assumed  to 
be  distributed  over  the  vibrational  energy  spectrum  with  a 
Boitzmann  distribution  characterized  by  a  vibrational  tempera¬ 
ture  Tv  _ .  The  relaxation  of  the  vibrational  energy  is  computed 
from  a  relation  of  the  form 


dUv.  UV.(T)  -  UVi(Tv.) 


where  is  an  empirical  relaxation  time.  If  desired,  any  dia¬ 
tomic  species  may,  instead,  be  considered  \.o  be  in  vibrational 
equilibrium  at  the  translational  temperature  T. 

Coupling  of  Vibrational  Energy  to  Rate  of  Dissociation 

With  the  assumption  that  dissociation  occurs  from  every 
vibrational  level  at  a  ^ate  proportional  to  the  population  at 
that  level  and  to  the  number  of  particles  having  a  translational 
kinetic  energy  of  sufficient  magnitude  to  cause  dissociation, 
the  rate  coefficient  for  the  dissociation  process  is  computed 
from  the  rate  coefficient  at  vibrational  equilibrium  (at  the 
same  temperature)  modified  by  a  factor  depending  upon  T,  Tv . , 
and  the  characteristic  vibrational  temperature  of  the  species. 


This  is  the  so-called  CVD  model. 
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Provision  for  More  Elaborate  Vibrational  Relaxation  Models 
The  contribution  of  each  chemical  reaction  ~o  the  produc 
tion  rate  and  removal  rate  of  the  related  species  is  individually 
computed,  and  the  net  production  is  determined  only  after  all 
reactions  have  been  treated.  In  this  way  not  only  are  possible 
losses  of  significant  figures  concentrated  at  one  point,  but  the 
net  molar  production  and  removal  rates  are  individually  available 
for  more  elaborate  vibrational  relaxation  models,  the  CVDV  model 

for  example. 

Choice  of  Chemical  Reactions 

A  large  number  of  chemical  reactions  (up  to  fifty)  may  be 
provided  with  the  structure  easily  adaptable  to  more  if  desired. 
The  computations  may  be  carried  out  with  any  combination  of 
reactions  as  well  as  with  none  at  all.  Not  only  may  air  chem¬ 
istry  be  treated,  but  also  hydrogen  combustion,  etc. 

Choice  of  Species 

The  program  is  so  structured  that  any  combination  of 
species  (initially  up  to  twenty)  may  be  considered.  Single 
species  may  be  studied  as  well  as  a  full  complement.  New 
species  may  be  added  without  program  modification  by  simply 
adding  a  standard  format  of  the  applicable  physical  data  for 


the  species. 
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Choice  of  Process 

Both  duct  flow  (prescribed  area  process)  and  the  pre¬ 
scribed  pressure  process  are  available,  and  the  computations 
may  proceed  from  one  to  the  other  in  a  single  problem. 

Provision  for  Wall  Friction  and  Heat  Transfer 
The  effects  of  friction  and  heat  transfer  are  accounted 
for  by  means  of  a  conventional  wall  friction  factor  and  a  heat 
flux  parameter.  The  heat  flux  may  be  treated  as  a  function  of 
stream  temperature,  if  desired,  and  may  thus  accommodate  thermal 
radiation  as  well  as  other  modes  of  heat  transfer. 

Modification  of  Analysis  to  Avoid  Iteration 
Usually,  the  only  differential  equation  needed  for  the 
thermodynamic  properties  is  the  momentum  equation,  which  for 
rrictionless  flows  is  simply  the  Euler  equation,  dp  =  -pVdV. 

When  pressure  is  prescribed,  the  velocity  derivative  immediately 
follows,  and  the  enthalpy  may  be  computed.  The  temperature  is 
then  computed  from  the  caloric  equation  of  state  by  a  time- 
consuming  trial-and-error  procedure.  Here,  however,  the  equa¬ 
tions  are  reformulated  to  yield  a  differential  equation  for  the 
temperature  so  that  the  enthalpy  may  now  be  determined  directly 
from  the  caloric  equation  of  state  without  iteration.  The  veloc¬ 
ity  is  then  found  from  an  algebraic  relationship. 
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Statistical  Calculation  of  Species  Data 

The  internal  energy  of  each  species  as  well  as  the  equi¬ 
librium  constants  for  all  reactions  are  calculated  from  parti¬ 
tion  functions  obtained  through  the  application  of  statistical 
mechanics  to  tabulated  spectroscopic  data. 

Flexibility  in  Treatment  of  Electronic  Energy 

The  electronic  energy  is  included  in  the  computet  ion  of 
the  internal  energy  of  the  species?  however,  the  number  of 
lev  ,1s  of  excitation  may  be  separately  chosen  for  the  individual 
species.  The  electronic  energy  may  also  be  omitted  entirely, 
if  desired.  Furthermore,  a  separate  electronic  temperature  is 
included,  common  to  all  species  as  a  first  approximation,  so 
that  the  electronic  energy  may  be  relaxed  if  a  suitable  relaxa¬ 
tion  differential  equation  is  provided. 

Provision  for  Oblique  Shocks 

An  oblique  shock  may  be  introduced  at  any  point  in  the 
calculation.  Either  the  shock  angle  or  the  turning  angle  may 
be  specified,  and  the  vibrational  energy  may  be  treated  as 
frozen  or  in  equilibrium  as  desired.  The  composition,  however, 
is  frozen. 

Comprehensive  Output 

The  output  of  data  as  the  calculations  proceed  is  designed 
to  provide  a  uprehensive  survey  of  not  only  the  gas  mixture, 
but  the  individual  species  and  the  chemical  reactions. 
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In  its  current  operational  form,  the  program  incorporates 
a  .39  reaction  scheme  for  13  pure  air  species  (0a ,  0,  Na  ,  N,  NO, 
e  ,  N0+,  03+,  0+,  Na+,  ,  Oa  ,  0  )  coupled  with  a  12  reaction 
scheme  for  the  species  Na,  Cl ,  NaCl,  NaO,  Na+,  and  Cl  .  The 
reactions  considered  and  their  associated  rate  constants  are 
given  in  Appendix  I . 


C - 3  AIR  CHEMISTRY  MODEL  WITH  THE  INCLUSION  OF  IMPURITIES 


The  one-dimensional  chemically  reacting  flow  program  dis¬ 
cussed  in  Section  C-l  has  been  used  to  simulate  the  experimental 
shock  tube  conditions  described  in  this  report.  The  calcula¬ 
tion  for  pure  air  under  typical  shock  conditions  (U  =  3  mm/jisec 

s 

Pi  =  2  mm  Hg;  T*  =  273°K),  utilizing  a  seven  species,  seven 
reaction  model  (0a ,  0,  Na ,  N,  NO,  e  ,  N0+,  reactions  1-7,  Appen¬ 
dix  I),  yields  the  electron  density  profile  shown  in  Fig.  C-3.1A 
A  nine  species,  nine  reaction  model  (add  to  above:  0a  ,  0  ; 
reactions  36,  37)  yields  essentially  identical  results. 

In  order  to  investigate  the  effects  of  some  naturally- 
occurring  contaminants  on  the  electron  density  profile  xn  other¬ 
wise  pure  air,  calculations  were  carried  out  with  trace  amounts 
of  NaC-L  and  Na.  The  initial  mass  fraction  of  contaminants  was 
varied  from  lO"*9  to  10-6 ,  and  the  resultant  electron  density 
profiles  arc  shown  in  Figs.  C-3.1B-, C  and  C-3.2.  Note  that  in 
all  cases  except  for  the  10”s  contaminant  mass  fraction,  where 
the  difference  is  slight,  the  equilibrium  electron  density 
level  is  equal  to  that  for  pure  air. 

From  Fig.  C-3.1  it  can  be  seen  that  for  pure  air,  or  for 
air  contaminated  with  NaC-t,  the  seven  species,  seven  reaction 
and  the  nine  species,  nine  reaction  air  models  yield  identical 
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results,  but  that  for  air  contaminated  with  Na  the  "nine  air" 
model  must  be  used.  This  is  a  consequence  of  the  varying  im¬ 
portance  of  different  electron-producing  reactions.  Figures 
C-3.3  and  C-3.4  depict  the  electron  production  rates  of  vari¬ 
ous  reactions  used  in  the  "nine  air"  plus  NaC^C.  and  "nine  specie  air" 
plus  Na  calculations.  It  can  be  seen  that  when  NaC>t  is  present, 
the  initial  production  of  electrons  behind  the  shock  is  due  to 
Ct  detachment,  and  reactions  involving  O3  and  0  are  nowhere 
of  importance.  On  the  other  hand,  when  Na  is  the  only  con¬ 
taminant,  Oa  detachment  accounts  for  the  initial  electron  pro¬ 
duction  mechanism.  In  both  situations,  far  downstream  of  the 
shock,  the  only  important  electron  producing  reaction  is  that 
of  N  fO  collisional  ionization. 
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C-4  SPECTROSCOPIC  SEARCH  FOR  IMPURITIES 

The  validity  of  an  analytical  modeL  used  to  represent  a 
physical  system  is  premised  on  the  fact  that  the  components  of 
the  physical  system  are  substantially  included  in  the  model. 

When  the  data  obtained  from  air  shock  is  compared  with  the 
model,  obviously  both  the  test  conditions  and  the  model  should 
represent  the  same  physical  system.  It  is  clear  from  the  re¬ 
cent  work  described  in  Section  C-3  that  small  amounts  of  impuri¬ 
ties  may  totally  mask  the  behavior  of  the  "pure"  system. 

To  investigate  what  impurities  do  exist  in  the  shock  tuoe 
system,  an  extensive  spectroscopic  study  was  undertaken.  The 
fundamental  tool  used  was  a  Jerrell  Ash  model  75,000  grating 
spectrograph.  This  instrument  has  an  extremely  large  effective 
aperture  (f~4).  The  grating  used  provided  approximately  20  &/mm 
dispersion  over  an  exit  width  of  about  15  cm. 

Extensive  calibration  of  the  instrument  has  been  made 
using  an  rf-excited  discharge  tube.  This  discharge  tube  was 
designed  so  that  any  given  gas  or  impurity  could  be  introduced 
and  tested  at  known  pressures.  While  the  conditions  of  this 
source  are  not  the  same  as  the  shock  tube,  the  basic  energy 
levels  of  the  species  are,  of  course,  unaltered,  even  if  the 
relative  intensities  are  dissimilar.  Hence  if  a  specific 
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impurity  is  suspected,  extensive  tests  can  be  made  to  establish 
its  detailed  spectra  which  can  then  be  used  to  verify  or  elimi¬ 
nate  the  existence  of  such  an  impurity  in  th  j  shocked  gas 
spectra. 

A  subordinate  piece  of  equipment  to  the  spectrograph  is  a 
recording  densitometer.  A  somewhat  unique  and  inexpensive 
densitometer  has  been  built  using  the  basic  parts  of  the  pre¬ 
cision  enlarging  facility  used  to  obtain  numerical  data  from 
the  scope  trace. 

The  spectrographic  plate  (5x7)  is  loaded  in  a  modified 
plate  holder  for  the  3x4  projector.  One-half  of  the  seven-, 
inch  plate  is  enlarged  about  15  diameters.  This  enlargement 
reduces  the  requirement  for  high  precision  stages  and  slits 
usually  required  in  such  an  instrument.  At  the  normal  image 
focal  plane,  a  set  of  silicon  photo  cells,  one  with  a  slit,  is 
mounted  on  a  precision  platform  driven  by  a  synchronous  motor. 

At  the  lens  of  the  projector,  a  light  chopper  is  used  to  achieve 
about  200  cps  fluctuation.  The  use  of  the  fluctuating  light 
source  allows  the  use  of  narrow  band  ac  a lif j ers  to  increase 
the  output  from  the  photo  cells,  thus  simplifying  the  design  of 
the  amplifiers  and  minimizing  the  effects  of  changes  in  ambient 
lighting  conditions. 
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Two  ohoto  cells  are  used,  one  with  a  narrow  slit  (.010")  to 
sample  the  trace,  ana  the  other  with  no  slit  (V  sq.  in  area)  to 
sample  the  unexposed  part  of  the  plate.  The  outputs  from  the 
amplifiers  for  those  two  photo  cells  are  put  in  bucking  series 
in  a  diode  detecting  circuit.  The  purpose  of  the  two  photo  cells 
is  to  compensate  for  changes  in  "fogging"  conditions  and  remove 
zero  drift  due  to  changes  in  light  intensity.  In  actual  use  no 
discernable  drift  has  been  observed  even  for  intensity  changes  of 
a  factor  of  two.  The  detected  compensated  output  is  recorded  or 
a  10"  synchronous  driven  strip  recorder. 

The  overall  system  from  spectrograph  to  output  chart  is  re¬ 
peatable  to  better  than  .1%  enabling  an  absolute  line  location  to 
better  than  3.0  A  out  of  3000  A,  which  is  twice  the  minimum  reso¬ 
lution  of  the  spectrograph. 

Three  specific  observations  of  the  emission  spectra  produced 
by  the  shocked  gas  were  made.  First,  an  "end  on"  measurement  was 
obtained  by  placing  a  window  at  the  end  of  the  rhock  tube.  Second, 
a  side  view  was  attempted  using  a  window  upstream  of  end.  Third, 
a  photoelectric  measurement  was  made  which  provided  the  time 
history  of  the  emission  behind  the  incident  shock. 

"End  on"Observation 

The  most  easily  obtained  integrated  spectra  is  provided  by 
terminating  the  shock  tube  with  a  window  and  collecting  the 
light  emitted  from  the  incident  and  reflected  shock.  The  in¬ 
tensity  of  the  light  thus  produced  is  sufficient  so  that 


the  event  may  be  easily  recorded  on  Tri  X  film  with  adequate 
resolution.  Unfortunately,  approximately  99%  of  the  recorded 
intensity  is  due  to  the  reflected  shock.  The  unfortunate  as¬ 
pect  is  presented  by  the  fact  that  the  emission  spectra  con¬ 
tains  species  that  may  have  been  introduced  from  the  walls  and 
window  by  mixing.  In  addition  the  interface  and  driven  gas 
may  also  be  mixed  and  heated  so  that  emission  may  occur  from 
the  gases  and  impurities  carried  down  the  tube. 

With  all  these  unsatisfactory  conditions,  the  overriding 
problem  of  adequate  light  is  satisfactorily  solved  and  the  in¬ 
formation  obtained  represents  the  total  aggregate  of  species  in 
the  system.  While  this  provides  broad  limits  on  the  total  im¬ 
purities,  no  correlation  can  be  obtained  as  to  if  the  species 
are  located  directly  behind  the  incident  shock,  the  region 
important  to  the  electron  density  profile  measurement. 

The  basic  experimental  arrangement  was  performed  to  mini¬ 
mize  some  of  the  emission  problem  discussed  above.  The  entrance 
lens  for  the  spectrograph  was  placed  so  that  its  focal  point 
was  at  the  inside  surface  of  the  end  window.  Secondly,  a  sys¬ 
tem  of  restricting  tubes  was  used  so  that  the  spectrograph 
admitted  light  principally  from  the  central  part  of  the  shock 


t  ibe. 
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The  impurities  observed  using  this  technique  include  Na, 

C,  Fe,  Ni,  Co,  Mn.  The  strong  presence  of  the  metallic  elements 
was  further  investigated  by  the  following  technique.  The  spec¬ 
tra  from  an  arc  using  stainless  steel  electrodes  was  recorded. 

A  negative  of  this  spectra  with  a  sodium  line  for  location  was 
printed.  This  "mask"  was  placed  in  the  exit  plane  of  the  spec¬ 
trograph  and  a  photo  cell  placed  behind  the  mask.  A  second 
photo  cell  was  placed  to  record  the  incident  light  to  the  spec¬ 
trograph.  An  electrical  difference  between  the  photo  cells  was 
formed.  After  suitable  gain  adjustments,  the  results  clearly 
indicate  that  primarily  all  of  the  stainless  steel  "emission" 
was  associated  with  the  arrival  of  the  shock  interface. 

To  avoid  the  mixing  problems  encountered  in  the  stagna¬ 
tion  region,  several  attempts  were  made  to  record  with  film  the 
spectra  emitted  by  the  plasma  directly  behind  the  incident  shock 
by  using  a  window  on  the  side  of  the  tube.  The  major  and  de¬ 
feating  problem  encountered  was  the  extremely  low  levels  of 
light  intensity.  Even  by  using  a  glass  tube  section  so  that 
light  could  be  gathered  from  the  largest  possible  area,  and  by 
processing  Tri  X  film  to  an  estimated  ASA  speed  of  2000,  no 
usable  image  could  be  obtained.  By  making  a  multiple  exposure 
of  thirty, shot  with  the  best  conditions,  only  faint  Na  lines 


were  visible. 
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As  a  final  solution  to  the  light  problem  the  spectrograph, 
was  instrumen-.  cd  with  5  sensitive  photomultipliers  (1P28)  .  The 
photomultipliers  were  positioned  behind  a  series  of  slits 
located  in  the  exit  plane  of  the  spectrograph.  The  slits  were 
made  adjustable  to  cover  the  equivalent  of  from  1  to  10  The 

housing  and  sockets  for  the  photomultipliers  were  provided  with 
circulation  passages  and  insulation  so  that  the  entire  ensemble 
could  be  cooled  toward  liquid  nitrogen  temperatures  to  reduce 
noise. 

Because  of  the  extreme  sensitivity  problems  anticipated, 
a  special  stabilized  power  supply  and  3  decade  logarithmic  am¬ 
plifiers  were  constructed.  The  somewhat  unique  duty  cycle 
requirements  of  this  particular  experiment,  permitted  i hf  de¬ 
sign  of  relatively  simple  logarithmic  amplifiers.  However,  the 
rise  time  requirement  of  better  than  3  jisec  severely  limited 
the  maximum  value  of  the  pnotomultipliers'  load  resistors.  A 
parametric  study  was  made  to  maximize  the  signal  (gain-bandwidth 
product)  co  the  noise.  The  amplifier  circuits  are  shown  in 
Fig.  C-4.1. 

This  wavelength  amplitude  calibration  was  accomplished  by 
using  a  standard  quartz  envelope  tungsten  libbon  lamp  and  a 
light  chopper.  After  computing  tie  emissivity  connected  black 
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body  radiation  curve,  the  individual  photomultipliers  were  checked 

over  the  spectral  range.  By  interchanging  individual  tubes  to 

compensate  for  their  particular  sensitivities,  an  arrangement  was 

» 

found  that  closely  matched  the  expected  black  body  curve  except 
for  the  red  region  as  shown  in  Fig.  C-4.2.  Because  of  the  in¬ 
herent  deficiency  of  the  1P28  photomultipliers  in  the  red  region 
(about  5900  A) ,  the  factor  of  ten  compromise  was  made. 

For  simplicity  and  interchangeability  of  amplifiers,  the 
gain  adjustment  was  ultimately  made  by  small  variations  in  the 
photomultipliers  load  resistors. 

The  wavelength  vs.  grating  position  for  each  slit  was 
achieved  using  a  low  pressure  Hg  discharge  tube.  This  light 
from  the  tube  was  chopped  (synchronously  with  the  power  line) 
and  the  output  from  the  photomultipliers  recorded  by  a  strip 
recorder  (one  channel  at  a  time) .  The  grating  was  slowly  ro¬ 
tated  with  a  gear- reduced  synchronous  motor.  This  procedure 
minimized  backlash  and  jitter  in  the  grating  rotation  system. 

The  location  of  13  principal  lines  was  noted.  These  locations 
were  least  square-fit  to  a  third  order  equation  (nine  fold 
redundancy)  and  an  accurate  set  of  grating  angle  vs.  wavelength 
tables  for  each  slit  computed. 

The  shock  tube  measurement  using  the  photoelectrical ly 
instrumented  spectrometer  was  made  from  the  side  wall  upstream 
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of  the  end.  The  window  was  contoured  quartz  about  2  cm  diam¬ 
eter.  It  was  found  necessary  to  liquid  nitrogen  cool  the  photo¬ 
multipliers  to  reduce  the  noise  as  anticipated.  The  shock  con¬ 
ditions  were- maintained  M~9,  1  mm  Hg  initial  pressure.  The 
repeatability  of  the  test  conditions  were  within  1%  in  shock 
velocity  by  using  the  mechanical  valve  as  driver.  An  auxiliary 
photo-sensitive  cell  was  used  to  monitor  the  total  light  output. 
The  integrated  power  recorded  ranged  typically  within  5%  of  the 
average . 

Because  of  the  multiplicity  of  lines,  it  was  necessary  to 
reduce  the  effective  exit  slit  width  to  2  On  this  basis,  to 

cover  the  complete  spectra  about  250  shots  would  have  been  re¬ 
quired.  Unfortunately  time  would  not  permit  the  complete  data 
recording  process  with  the  accompanying  enormous  data  reduction 
that  would  be  required  (analysis  of  1250  pictures) .  Instead, 
selective  regions  of  the  spectra  were  chosen  for  detailed  analy¬ 
sis.  Some  of  these  selected  lines  ^2  %  wide)  are  shown  in 


Fig.  C-4.3. 


Fig.  C-4.3  Selected  Emission  Lines  Observed 
Behind  Incident  Shock 


Each  curve  represents  the  mean  behavior 
taken  for  the  two  lines  indicated  for  an 
average  of  two  shots  each. 
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C-5  COMMENTS  ON  THE  EXPERIMENTAL  AND  COMPUTED  ELECTRON 

DENSITY  PROFILE 

The  basic  shock  conditions  for  which  this  comparison  has 
been  made  were  detailed  in  Section  B-3.  The  specific  velocity  and 
initial  pressure  correspond  to  Mach  8.9  shock  (3025  m/s)  at  ini¬ 
tial  pressure  of  134,000  feet  altitude  (2.0  mm  Hg) .  The  experi¬ 
mental  curve  in  Fig.  C-5.1  was  taken  from  the  "mean"  determina¬ 
tion  discussed  in  Section  B-3. 3. 

The  analytically  determined  curve  in  Fig.  C-5.1  was  cal¬ 
culated  by  the  program  discussed  in  Section  C-l.  The  obvious 
order  of  magnitude  difference  in  Fig.  C-5.1  warrants  serious 
consideration . 

The  error  indicated  either  lies  with  the  analytical  model 
used  for  the  computed  curve  or  with  the  experimental  technique 
used  to  measure  the  electron  density  profile,  or  of  course, 
some  combination  of  each. 

The  experimental  technique  is  discussed  in  detail  in  Sec¬ 
tion  B  and  in  particular.  Section  B-4  provides  such  support¬ 
ing  evidence  as  to  validate  the  result  within  the  broad  con¬ 
text  of  the  following  discussion. 

The  analytical  model  may  be  incorrect  because  (1)  the 
constituents  of  the  model  and  that  of  the  experiment  differ. 
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(2)  seme  rate  or  energy  data  is  incorrect  or  (3)  the  computa¬ 
tional  procedure  is  invalid.  It  is  probable  that  the  error  is 
attributable  to  (1)  even  though  (2)  or  ’,3)  cannot  be  totally 
eliminated.  In  reference  to  item  (2)  several  parametric  studies 
have  been  made  to  ascertain  the  sensitivities  of  the  various 
reactions  and  reacting  systems.  These  studies  showed  that  no 
reasonable  alteration  of  constants  could  be  made  to  satisfy 
the  results  of  both  the  local  test  conditions  of  this  experi¬ 
ment  and  other  data  for  higher  temperature  conditions  where  the 
numerical  analysis  is  in  much  better  agreement  with  experiment. 

On  the  other  hand  in  reference  to  item  -1)  above,  the  ex¬ 
clusion  of  an  impurities  reacting  system  could  reconcile  the 
apparent  difference.  In  Section  C-2,  the  influence  of  sodium 
on  the  electron  rate  production  is  shown.  It  is  clear  from  the 
curves  in  Fig.  C-3.2  that  for  concentrations  of  salt  similar  to 
that  approximating  natural  occurrence  near  sea  water  (1  part  per 
10s).  the  electron  profile  should  be  extremely  steep  and  yet 
it  may  be  noted  that  the  equilibrium  levels  are  essentially  un¬ 
altered.  Since  the  experimental  results  easily  fall  between 
the  "pure"  air  case  and  the  slightly  contaminated  cases,  it  is 
likely  that  sodium  or  a  similar  agent  could  account  for  the  dif¬ 


ference  observed. 
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It  is  worthy  to  note  that  the  results  presented  in  Sec¬ 
tion  C-2  assume  that  the  sodium  occurs  uniformly  in  the  gas 
prior  to  the  shock  arrival.  From  the  time  resolved  spectro¬ 
scopic  data  in  Fig.  C-4.3A  the  initial  threshold  observance  of 
the  sodium  emission  occurs  15  cm  behind  the  shock  front  and  then 
increases  rapidly.  In  addition  it  may  be  noted  that  the  iron 
emission  follows  roughly  the  same  behavior  at  a  later  time,  but 
before  the  shock  interface.  The  suggestion  of  this  data  is 
that  the  sodium, like  the  iron  , is  coming  from  the  walls  and  is 
mixed  and  diffused  into  the  shocked  region.  With  this  mechan¬ 
ism  working,  the  initial  reacting  region  behind  the  shock  could 
be  behaving  as  the  "pure  air"  model  suggests  while  the  region 
slightly  upstream  of  the  shock  frcnt  may  be  influenced  by  the 
sodium  diffusing  and  mixing  from  the  walls.  One  piece  of  evi¬ 
dence  missing  from  this  analysis  is  the  threshold  sensitivity  of 
the  spectrographic  system  to  sodium.  Unfortunately  the  experi¬ 
ment  evaluation  of  this  threshold  is  complicated  and  could  not 
be  attempted. 
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CHAPTER  V 

THEORETICAL  STUDIES  ON  NON-UNIFORM  PLASMAS 

D-l  INTRODUCTION 

The  following  chapter  contains  a  summary  and  review  of 
some  of  the  analytical  work  generated  as  an  adjunct  to  the  ex¬ 
perimental  work.  The  first  section  considers  the  general  analy¬ 
sis  of  the  electromagnetic  properties  of  non-uniform  plasmas 
from  microscopic  and  macroscopic  considerations.  Section  two 
presents  a  detailed  analysis  of  diffusion  effects  in  a  non- 
uniform  plasma.  Section  three  considers  the  electromagnetic 
propagation  characteristics  of  a  shock  tube  produced  plasma  with 
a  radially  non-uniform  plasma  density. 

It  is  well  known  that  in  a  first  order  approximation,,  the 
propagation  of  an  electromagnetic  wave  may  be  analyzed  with  a 
simplified  scheme  in  which  the  local  properties  of  the  ionized 
gas  are  given  by  the  plasma  frequency  and  the  constant  colli¬ 
sion  frequency.  Even  in  this  case  the  solution  of  two-  or 
three-dimensional  propagation  problems  may  present  tremendous 
mathematical  difficulties.  On  the  other  hand,  in  a  general 
case  for  a  non-uniform  gas,  the  effort  required  to  solve  these 
problems  is  not  justified  due  to  the  oversimplified  scheme  of 
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the  electrical  properties  of  the  gas.  When  the  change  of  the 
macroscopic  thermodynamic  properties  of  the  gas  becomes  important 
in  a  length  of  the  same  order  of  magnitude  of  a  wavelength  of 
the  electromagnetic  field,  it  becomes  necessary  to  analyze  the 
possible  effects  of  additional  phenomena  such  as  the  transport 
process,  the  anisotropic  behavior  of  the  electrical  properties, 
the  electrical  polarization  of  the  medium,  etc.  Furthermore, 
the  coupling  between  transverse  oscillations  and  coherent  longi¬ 
tudinal  electron  oscillations  may  become  a  significant  factor 
leading  to  a  substantial  increase  of  the  dissipation  of  electro¬ 
magnetic  energy  in  thermal  energy  in  the  ionized  gas.  If  the 
relative  importance  of  these  phenomena  could  be  properly  eval¬ 
uated  from  the  theoretical  standpoint,  then  it  would  be  possible 
to  define  the  set  of  parameters  which  govern  the  behavior  of  an 
ionized  gas  for  each  particular  case.  The  major  difficulty  of 
a  theoretical  approach  is  due  to  the  fact  that  it  requires  de¬ 
tailed  information  about  the  microscopic  processes,  whic'  at 
the  present  time  is  rather  incomplete,  particularly  for  the 
case  of  weakly  ionized  air  as  encountered  in  flight  conditions. 
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D-2  ELECTROMAGNETIC  PROPERTIES  OF  NONUNIFORM  PLASMAS 
{ 1 )  Microscopic  Considerations  Ref.  4) 

The  effect  of  microscopic  inhomogeneities  on  the  elec¬ 
tromagnetic  characteristics  of  a  plasma  are  of  interest.  In 
order  to  estimate  such  effects  it.  is  necessary  to  calculate 
curtain  distribution  functions  directly  deducible  from 
Boltzmann's  equations.  If  inhomogeneities  are  present  (due  to 
c  "adients  of  concentration,  of  temperature,  etc.)  as  far  as 
we  know,  no  actually  developed  method  of  calculation  exists. 

For  this  reason  we  decided  to  examine  various  models 
beginning  with  the  simplest  ones: 

(a)  Plasma  assumed  to  be  a  weakly  ionized 
Lorentzian  gas.  considering  electrons  and 
neutral  molecules  and  neglecting  Coulomb 
interactions . 

(b)  Three -component  plasma  -'electro  1 ,  ions,  neutral 
molecules)  taking  Coulomb  interactions  into 
account . 

Since  the  analysis  of  model  2  has  appeared  particu¬ 
larly  difficult  from  a  general  point,  of  v; ew  (for  the  reason  it 
implies  the  solution  of  a  system  of  Boltzmann's  equations), 
first  of  all  we  decided  to  carry  out  in  a  detailed  way,  a 


method  of  actual  calculation,  in  order  to  face  the  simplest 
model  1  even  if  it  is  not  very  realistic.  In  addition,  for  the 
scheme  1  we  have  decided  to  examine  separately  the  various  in¬ 
homogeneity  effects  and  also,  to  keep  near  to  the  model  studied 
from  the  macroscopic  viewpoint,  we  have  selected  the  detailed 
examination  of  the  influence  of  temperature  gradient  in  a  slab 
containing  the  plasma,  whose  walls  are  absorbing  and  emitting 
with  a  Maxwellian  distribution. 

The  complete  physico-mathematical  characteristics 
(in  our  opinion  not  so  wholly  developed  in  the  literature)  of 
model  1  are  described  up  to  the  deduction  of  final  equations. 

In  addition,  a  known  analytical  method  used  for  solv 
ing  the  equations  of  Section  1  (series  expansion  of  eigen  func 
tions  of  the  collision  operator)  is  described. 

(2)  Macroscopic  Properties  (Refs.  3  and  4) 

The  electromagnetic  properties  of  a  macroscopically 
neutral  non-uniform  plasma  are  considered  by  an  approach  based 
on  the  thermodynamics  of  irreversible  processes.  A  three-fluid 
model  of  the  plasma  is  assumed  end,  in  this  first  stage,  the 
plasma  is  considered  to  be  at  rest  in  a  suitable  reference 
frame,  the  components  are  in  mutual  thermal  equilibrium  and 
imposed  m  jnetic  fields  are  absent.  The  non-uniformities  are 
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those  connected  with  the  presence  of  gradients  of  state  param¬ 
eters  such  as  temperature,-  pressure  and  concentrations. 

The  presence  of  these  gradients  induces  a  polariza¬ 
tion  in  the  medium  even  in  the  absence  of  external  electric 
fields.  It  will  be  shown  that,  in  the  stationary  state,  the 
induced  polarization  depends  on  only  two  independent  gradients, 
which,  for  convenience,  are  taken  to  be  those  of  the  tempera¬ 
ture  and  the  pressure.  An  expression  for  the  induced  polariza¬ 
tion  is  derived,  which  is  valid  for  a  most  general  case  which 
is  later  simplified  to  the  case  of  the  imperfect  Lorentz  gas 
(i.e.,  weakly  ionized  plasma  for  which  the  macs  ratio  |  between 
the  negative  charge  carriers  and  neutral  molecules  is  much 
smaller  than  one)  and  the  perfect  Lorentz  gas  (i.e.  weakly 
ionized  plasma  with  £  =  0) . 

The  phenomenological  coefficients  relating  the  in¬ 
duced  polarization  to  the  temperature  and  pressure  gradients 
are  first  expressed  in  terms  of  thermodynamic  properties  of 
the  plasma.  It  is  then  shown  how  these  coefficients  can  be 
expressed  only  in  terms  of  simple  binary  diffusion  coefficients 
and  thermal  diffusion  coefficients  pertinent  tc  the  plasma  con¬ 
stituents.  This  last  step  furnishes  the  linking  element  be¬ 
tween  the  present  macroscopic  approach  and  the  mere  refined, 
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but  usually  less  general,  statistical  approach.  Indeed  it  indi¬ 
cates  how  the  subject  phenomenological  coefficients  can  be 
rigorously  evaluated,  in  terms  of  the  "microscopic"  character¬ 
istics  of  the  plasma  constituents,  once  the  statistical  approach 
has  gone  far  enough  to  provide  workable  results.  In  the  mean¬ 
time,  nowever;  the  fact  that  one  needs  to  know  only  the  binary 
transport  coefficients  provides  a  possibility  of  readily  per¬ 
forming  order  of  magnitude  analysis  by  utilizing  the  available 
experimental  data  on  electron-ions,  electron-molecules  and  ion- 
molecules  collisions. 

The  results  of  the  prese,  *-  stage  of  the  analysis 
already  open  up  two  possible  ways  of  utilization,  both  of 
noticeable  interest. 

In  the  first  place,  as  mentioned  before,  several 
order-of-magnitude  comparative  analyses  can  be  carried  out. 

For  instance,  one  can  evaluate  the  relative  importance  of  the 
different  inhomogeneities  with  respect  to  the  polarizations 
they  are  able  to  induce  (e.g..  is  a  temperature  or  a  pressure 
gradient  more  important  to  this  effect?) .  Also,  one  can  compare 
the  polarization  induced  by  the  inhomogeneities  with  that  in¬ 
duced  by  an  external  electric  field  and  thus  determine  the 
ranges  within  which  the  former  one  can  be  neglected  and  the 
medium  treated  as  homogeneous. 
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The  second  application  of  equal  practical  interest 
follows  from  the  consideration  of  the  electromagnetic  field 
which  is  created  by  the  induced  polarization  (which,  in  princi¬ 
ple,  can  be  computed  by  evaluating  the  Hertz  potential  whose 
source  is  the  induced  polarization) .  When  we  induce  into  the 
plasma  a  polarization  by  means  of  controlled  suitable  inhomo¬ 
geneities,  part  of  the  energy  furnished  to  maintain  the  gradi¬ 
ents  is  spent  to  "create"  the  induced  electromagnetic  field. 
This  is  nothing  but  the  underlying  idea  of  energy  conversion 
systems  and  the  present  results  can  be  used  to  investigate  both 
the  feasibility  and  efficiency  of  a  number  of  such  systems. 

The  steps  of  the  analysis  reported  are  as  follows: 

The  irreversible  thermodynamic  description  of  the 
system  is  performed.  The  extensive  and  intensive  state  param¬ 
eters  and  the  pertinent  mass  and  energy  fluxes  for  a  mixture  of 
three  fluids,  of  which  two  have  negative  and  positive  charges, 
are  defined.  The  "kinetic"  relations  between  fluxes  and  gen¬ 
eralized  forces  are  established  and,  through  suitable  use  of 
the  basic  theorems  of  the  thermodynamics  of  irreversible  proc¬ 
esses,  the  general  expression  for  the  mass  fluxes  in  terms  of 
pressure;,  temperature  the  electric  potential  ($)  gradients  and 
thermodynamic  properties  of  the  plasma  is  arrived  at. 
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The  general  expression  for  the  polarization  is  ob¬ 
tained,  for  a  macroscopically  neutral  plasma,  in  the  form: 

grad  $  =  grad  T  +  grad  p 

where  the  cr  are  phenomenological  coefficients.  Their  expres¬ 
sion  in  terms  of  the  thermodynamic  properties  of  the  sma  is 

also  derived  in  this  section,  both  for  the  general  case  and  for 
the  simplified  case  of  the  Lorentz  gas. 

The  problem  of  relating  the  coefficients  or  to  binary 
transport  coefficient  is  considered.  By  suitable  transformation 
of  fluxes  and  affinities  it  is  shown  how  they  can  be  expressed 

in  terms  of  three  binary  diffusion  coefficients  Dia,  Da3,  D13 

T  T 

and  two  thermal  diffusion  coefficients,  Dx * ,  Da  ,  which  refer 

to  the  plasma  constituents.  General  expressions  are  presented 

but  no  attempt  is  yet  made  to  either  simplify  them  or  evaluate 

their  comparative  order  of  magnitude. 

(3)  Electromagnetic  Properties  of  Non-uniform  Plasmas 
in  Thermal  Equilibrium,  Isotropic  Case  (Ref.  3) 

Using  the  three  fluid  models  and  an  approach  based 
on  the  thermodynamics  of  irreversible  processes,  the  effect  of 
non-uniformities  due  to  gradients  of  state  parameters  (tempera¬ 
ture  and  density)  on  the  electromagnetic  properties  are  consid¬ 
ered.  The  presence  of  these  gradien  s  induces  an  electric  fi^ld 
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(i.e.f.)  in  the  medium,  even  in  the  absence  of  external  electric 
fields.  It  will  be  shown  that,  in  the  stationary  state  and  in 
conditions  of  macroscopic  neutrality,  the  induced  e.f.  depends 
only  on  one  independent  gradient  which,  for  convenience,  is 
taken  to  be  that  of  temperature.  An  expression  for  the  i.e.f. 
is  derived,  which  is  valid  for  a  most  general  medium  and  which 
is  later  simplified  to  the  case  of  an  imperfect  Lorentz  gas 
(i.e.,  weakly  ionized  plasma  for  which  the  mass  ratio  between 
the  negative  charge  and  neutral  molecules  is  much  smaller  than 
one)  . 

The  electrical  conductivity  <j  defined  as  the  ratio 
between  the  electric  current  I  and  the  electric  potential  gradi¬ 
ent  in  a  first-order  stationary  state  in  which  a  given  non- 
uniform  distribution  of  V  $  is  maintained  on  the  boundaries  of 
a  system  is  also  determined. 

The  phenomenological  coefficients  are  first  expressed 
in  terms  of  thermodynamic  properties  of  the  plasma.  It  is 
shown  how  these  coefficients  can  be  expressed  in  terms  of  only 
simple  binary  molecular  diffusion  coefficients  and  thermal  dif¬ 
fusion  coefficients  pertinent  to  the  plasma  constituents. 

This  last  step  makes  it  possible,  by  utilizing  the 
available  experimental  and/or  theoretical  data  on  electron- ion. 
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electron-molecule,  and  ion-molecule  collisions,  to  express  the 
phenomenological  coefficients  in  terms  of  the  local  thermo¬ 
dynamic  state  of  the  medium.  It  was  possible  to  observe  that 
the  i.e.f.  increased  with  the  temperature  difference  and  the 
degree  of  ionization.  These  results  qualitatively  agree  with 
the  theoretical  conclusions  which  are  arrived  at  in  the  present 
report. 

The  steps  of  the  analysis  are  as  follows:  Part  1, 
the  irreversible  thermodynamic  description  of  the  system  is  per¬ 
formed.  The  extensive  and  intensive  state  parameters  and  the 
pertinent  mass  and  energy  fluxes  for  a  mixture  of  three  fluids, 
one  with  negative  and  one  with  positive  charges,  are  defined. 

The  dynamic  relations  between  fluxes  and  generalized  forces  are 
established  and,  through  suitable  use  of  the  basic  theorems  of 
the  thermodynamics  of  irreversible  processes,  the  general  ex¬ 
pression  for  the  mass  fluxes  in  terms  of  electron  concentration, 
temperature,  and  electric  potential  gradients  is  arrived  at. 
Parts  2  and  4,  the  general  expressions  for  the  i.e.f.  and  the 
electrical  conductivity  are  obtained.  Part  3,  the  problem  of 
relating  these  phenomenological  coefficients  to  binary  transport 
coefficients  is  considered.  By  suitable  transformation  of 
fluxes  anC  affinities  it  is  shown  how  they  can  be  expressed 
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in  terms  of  three  binary  diffusion  coefficients  Dia ,  Das#  D13 

T  T 

and  two  thermal  diffusion  coefficients  Dx  ,  Ds  ,  which  refer 
to  the  plasma  constituents. 

In  addition,  the  evaluation  of  the  molecular  and 
thermal  diffusion  coefficients  for  a  plasma  is  carried  out  and 
in  Part  2  an  analysis  of  the  order  of  magnitude  for  the  trans¬ 
port  coefficients  is  performed  for  the  case  of  a  weakly  ionized 
imperfect  Lorentz's  gas. 
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D-3  DIFFUSION  EFFECTS  IN  A  NON-UNIFORM  PLASMA 
(1)  Basic  Equations 

In  the  experimental  conditions  encountered  in  the 
shock  tube  measurement  the  gas  is  still  weakly  ionized.  There¬ 
fore,  in  a  theory  describing  the  electric  properties  of  the  gas, 
the  collisions  between  charged  particles  may  be  neglected  and 
only  electron-neutral  collisions  and  ion-neutral  collisions 
will  be  taken  into  account.  Furthermore,  one  neglects  the  per¬ 
turbation  of  the  neutral  component  of  the  gas  due  to  the  ioniza¬ 
tion  and  recombination  processes  and  to  the  motion  of  the  charged 
particles.  Thus  the  electron  and  ion  momentum  equations  are 
written  in  the  form: 
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In  Eq.  (D-l)  no  external  electromagnetic  field  is  ap¬ 
plied  to  the  plasma  and  E  is  the  intensity  of  electric  field 

—4  ^ 

induced  by  a  small  space  charge  distribution;  u.f  u  ,  u  ,  are 

i  e  n 

the  ion,  electron  and  neutral  particles  average  velocities;  n^, 

n  ,  are  the  ion  and  electron  densities  and  q  is  the  value  of 
e 

the  elementary  electric  charae,  m,,  m  ,  m  ,  are  the  ion,  elec- 

l  e  n 

tron  and  neutral  particle  masses  and 
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V is  the  collision  frequency  of  an  electron  with  the  neutral 
particles  and  v ^  is  the  collision  frequency  of  an  ion  with  the 
neutral  particles.  One  assumes  that  electrons  and  ions  are  in 
thermal  equilibrium  with  the  neutral  gas.  Thus  if  the  tempera¬ 
ture  T  is  assumed  constant  throughout  the  plasma,  one  has: 


vp  =  k  ?  V  n  ,  V  p,  -  k  T  v  n. 
e  e  i  i 


(D-3) 


where  k  is  the  Boltzmann  constant.  The  left  side  of  Eq.  (D-l) 
corresponds  to  the  electron  and  ion  accelerations.  Let  us 
confine  the  analysis  to  a  steady  state  situation  where  d/dt  =  0. 
Outside  of  the  plasma  sheet  at  the  wall  of  the  container  and 
in  a  region  of  sufficiently  small  gradients  of  electron  and  ion 
densities,  the  acceleration  terms  are  negligible  compared  to  the 
last  terms  on  the  right  side  of  Eq.  <p-l) ,  and  the  momentum  equa¬ 
tions  may  be  written  again  in  the  form: 
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u  ,  U./  are  the  electron  and  ion  mobilities.  D  ,  D.  are  the 
^e  *i  e  l 


electron  and  ion  diffusion  coefficients.  The  intensity  of  elec¬ 


tric  field  is  related  to  the  electron  and  ion  densities  through 


the  Poisson's  equation 


7-E  =  (n.-n  )  (D-6 ) 

C  i  e 
o 


where  cQ  is  the  dielectric  permeability  or  permittivity  of 
vacuum.  Furthermore,  in  a  steady  state  situation  the  electron 
and  ion  continuity  equations  reduce  to 


V  (n  u  )  =  7*  (n  .  u  •  )  =  h  (D-7) 

e  e  ii  e 

where  n  indicates  the  rate  of  creation  of  electrons  per  unit 
e 

volume;  depends  upon  the  electrons,  ions  and  neutral  par¬ 
ticle  densities  according  to  the  particular  ionization  and 
recombination  processes  which  occur  in  the  ionized  gas.  By 
virtue  of  (D-7),  from  Eq.  (D-4)  one  obtains: 
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Outside  of  the  plasma  sheet  and  in  a  region  of  sufficiently 
small  density  gradient,  one  can  assume 


n  -n .  «  n 

e  l  e 


(D-9) 


Thus  the  sum  of  Eq.  (D-8) leads  to 
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where 
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D  is  the  amb-ipolar  diffusion  coefficient. 

The  local  electric  current  density  is  given  by 


j  =  q  (n.u. -  n  u  ) 
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By  virtue  of  Eqs.  (D-4)  and  (D-9),  Eq.  (D-12)  becomes: 
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E  may  be  written  in  terms  of  the  electrostatic  potential  $ 

E  =  -V$  (D-14 } 


In  particular  if  no  electric  current  is  flowing  in  the  plasma, 
and  no  external  electric  field  is  applied  to  the  gas,  Eq.  (D-I3) 
reduces  to 
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where  <f>Q  is  an  arbitrary  constant. 

From  Eq.  (D  13)  one  obtains 
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Hence  the  difference  between  electron  and  ion  densities  is: 
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It  is  of  interest  to  observe  that  Eq. (D-15)  implies 
that  in  the  absence  of  an  electric  current,  the  electrons  are 
trapped  in  the  distribution  of  electric  field  due  to  the  space 
charge  in  the  plasma. 

( 2 )  Boundary  Condiuions 

Inside  the  plasma  we  assume  that  the  electrostatic 
potential  is  zero  where  the  electron  density  becomes  equal  to 
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the  equilibrium  density  n_  .  is  the  electron  equilibrium 

eo  eo 

density  which  corresponds  to  the  uniform  temperature  and  density 
of  the  neutral  particles  of  the  gas.  Eq.  (D-15* becomes 

D-D.  n 

$  =  — — -  In  —  (D-18) 
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By  virtue  of  Eqs.  <D-5),  D  »  D.,  4  »  u.  .  Thus(0-181  reduces  to 

e  1  e  1 

ne  *  nr.  e  kT  (D-19 ) 
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Let  the  plasma  sheet  be  defined  as  the  small  layer 
adjacent  to  the  wall,  where  collisions  of  the  charged  particles 
with  the  neutral  particles  may  be  neglected.  On  the  other  hand 
a  strong  negative  potential  of  the  wall  with  respect  to  the 
plasma  must  arise  in  order  to  maintain  the  condition  of  zero 
electric  current.  Let  x  be  the  thickness  of  the  plasma  sheet, 
and  let  us  analyze  the  electron  and  ion  motion  in  the  steady 
state.  One  may  analyze  the  motion  assuming  a  one-dimensional 
problem  where  the  axis  x  is  perpendicular  to  the  wail  and  x=  0 
at  the  wall.  If  the  electrons  are  almost  entirely  reflected 
back  into  the  plasma  due  to  the  negative  potential  distribution 
in  the  proximity  of  the  wall,  throughout  the  plasma  sheet  solu¬ 
tion  (D-19)  is  valid,  as  far  as  the  electron  density  is  concerned. 
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On  the  other  hand  strong  accelerations  of  the  ions 
may  be  expected  in  the  plasma  sheet  and  the  ion  momentum  equa¬ 
tion  leads  to 
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where  u.  and  $  are  the  ion  velocity  and  the  electrostatic 
ixs  s 

potential  at  the  edge  of  the  plasma  sheet  (x=xs).  The  con¬ 
tinuity  of  ion  velocity  at  x =  xg  requires  that 
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where  (d$/dx)  is  the  derivative  of  $  at  x=x  .  Let  4>  be  the 
s  so 


electrostatic  potential  of  the  wall  and  assume  that  ions  and 
electrons  recombine  as  they  collide  with  the  wall.  The  number 
of  electrons  reaching  the  wall  per  unit  time  and  unit  area  is 
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On  the  other  hand  the  ion  flow  toward  the  wall  per  unit  area  and 
unit  time  is 


-  n  ,  u  ,  =  -  n , 

i  ix  is 


u , 
ixs 


qDn^ 

o 


kT 


dj$> 

dx 


(D-23) 


'  >s 

With  no  electric  current  at  the  wall  the  ion  flow  must  be  equal 
to  the  electron  flow.  Thus: 
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and  the  ion  velocity  at  x  =  x  is: 
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Thus  from  Eq.  (D-20)  one  obtains 
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Inside  the  plasma  sheet  the  Poisson  equation  reduces  to 
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which  may  be  easily  integrated,  leading  to 
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where  C  is  the  integration  constant.  The  constant  C  is  computed 
with  the  aid  of  the  boundary  condition 
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Without  entering  into  the  details  of  the  solution  of  Eq.  (D-29) 
one  may  observe  that  d2$/dxE  has  a  negative  or  zero  value  if 
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In  particular,  condition  (D-32)  is  verified  in  the  proximity  of 
the  edge  of  the  plasma  sheet  ( <£«  $  )  if 
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m . 

l 


2  rr  m 

e 


'IS 

jsT 


*  1 


(D-34) 


By  virtue  of  (D-25) ,  condition  (D-34)  becomes 
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If  condition  (D-35)  were  not  satisfied,  the  ion  density 
would  drop  faster  than  the  electron  density  at  the  entrance  of 
the  plasma  sheet  and  this  would  be  inconsistent  with  a  distribu¬ 
tion  of  electrostatic  potential  inside  the  plasma  sheet  where  $ 
steadily  decreases  toward  the  wall. 

( 3 )  Steady  State  Distribution  of  a  Plasma  in  a  Cylinder 

It  is  of  interest  to  discuss  in  some  detail  i he  c-u  . se¬ 
quences  of  the  previous  considerations  ::  a  situation  clou.  ~ar  to 
the  experimental  conditions  in  the  shock  tube.  In  order  to  do 
so,  assume  a  simple  case  for  which  the  recombina*  ’on  rate  of  the? 
electric  charges  is  negligible  and  the  electron  creation  rate 
may  be  written  in  the  form 

ne  =  a  ne  4D-36) 

where  a  is  a  function  of  the  neutral  particle  density  alone  for 
a  given  value  of  the  electron  temperature.  Equation  (D—  36 ) 
would  correspond  to  an  ionization  process  due  to  electron- 
neutral  particle  collisions.  Assume  further  a  cylindrical  sys¬ 
tem  of  reference  r,  <£,  z,  where  r  =  r  ,  corresponds  to  the  wall 

o 

of  the  cylinder.  Furthermore  the  neutral  particle  velocity  u^ 
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is  assumed  to  be  uniform  and  parallel  to  the  z-axis  and  u.  is 
assumed  to  be  constant.  From  Eq.  (D-10)  one  obtains 
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where  no  $  dependence  of  the  electron  distribution  is  taken  into 
account.  A  particular  solution  of  (D-37)  in  the  plasma  outside 
of  the  plasma  sheet  is: 
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where  z  satisfies  the  equation 
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J  is  the  Bessel  function  of  the  first  kind  and  order 
o 

zero  and  y  is  a  constant  to  be  determined  by  the  boundary  condi¬ 
tions.  Assume  that  a  steady  state  situation  is  found  where  the 
electron  density  is  independent  of  z.  In  this  case  the  constant 
y  is  referred  to  the  coefficient  a  and  the  ambipolar  diffusion 
coefficient  through  the  equation 
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(D-40) 
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The  radial  ion  velocity  is  given  by: 
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Let  rg  be  the  radius  of  the  edge  of  the  plasma  sheet.  At  the 

edge  of  the  plasma  sheet  the  space  charge  separation  is  negligible. 

Consequently#  by  virtue  of  condition  (D-35) #  the  minimum  value 

of  u,  at  r  =  r  must  be: 
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and  by  virtue  of  (D-41)  one  obtains: 
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Taking  into  account  the  definition  of  the  ambipolar 

diffusion  coefficient  and  assuming  that  the  thickness  r  -  r  of 

o  s 

the  plasma  sheet  is  small  compared  to  the  radius  of  the  cylinder# 
condition  (D-43)  leads  to 
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i .  e. 


Ji  (yrQ) 
yro  J0 (yrQ) 


where 


2_  I  2kT 

"i  J  mi 


(D-45) 


(D-46) 


is  the  mean  free  path  for  the  ions.  Taking  into  account  that 
must  be  small  compared  to  the  radius  of  the  cylinder,  Eq.  (D-45) 
reduces  to 

X. 

J0(yrQ)  *  2.01  ~  (D-47 ) 

o 

ne  must  be  a  positive  density  throughout  the  plasma.  Thus 

yr  <  2.405  and  with  X.  «  r  ,  Eq.  (D-47)  leads  to  the  conclu- 
o  l  o 

sion  that  the  electron  density  at  the  edge  of  the  plasma  sheet 
must  be  small  compared  to  the  electron  density  at  the  center  of 
the  cylinder.  It  is  worthwhile  pointing  out  that  condition 
(D-47)  implies  that  the  rate  of  creation  of  electrons  in  the 
plasma  balances  the  rate  at  which  electrons  and  ions  diffuse 
toward  the  wall  of  the  cylinder. 

(4)  Concluding  Remarks 

The  quantitative  solution  of  the  system  of  equations 
of  Section  (1)  and  of  Eq.  (D-29)  must  be  carried  out  in  order 
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to  determine  the  electron  density  inside  the  plasma  and  the 
thickness  of  the  plasma  sheet.  Nevertheless,  the  above  consid¬ 
erations  indicate  that  a  large  effect  in  the  transversal  diffu¬ 
sion  should  be  expected  due  to  the  large  value  of  the  diffusion 
velocity  of  the  ions  at  the  edge  of  the  plasma  sheet.  As  a  rela¬ 
tively  simple  illustration  of  the  above  phenomena,  the  one¬ 
dimensional  diffusion  effects  between  two  regions  of  air  at 
different  temperatures  and  separated  by  a  plane  boundary  are 
presented  in  Appendix  II. 
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D-4  ELECTROMAGNETIC  PROPAGATION  IN  A  WAVEGUIDE  CONTAINING 

A  RADIALLY  NON-UNIFORM  PLASMA 

The  object  cf  this  section  is  to  study  a  particular  prob¬ 
lem  of  propagation  of  electromagnetic  waves  in  a  plasma  contained 
inside  a  metallic  cylinder.  The  plasma  is  assumed  to  be  non- 
uniform  in  the  radial  direction,  but  its  properties  are  inde¬ 
pendent  of  both  the  axial  and  angular  coordinates. 

Specifically,  the  non-uniform  properties  of  the  plasma  cor¬ 
respond  to  an  electron  density  which  is  a  function  of  the  radial 
coordinate.  Furthermore,  we  assume  that  the  plasma  collision 
frequency  is  negligible  with  respect  to  the  frequency  of  the 
electromagnetic  wave.  Such  a  radial  non-uniformity  of  the 
plasma  can,  for  example,  be  attributed  to  the  ambipolar  diffu¬ 
sion  process  of  charges  which  move  perpendicular  to  the  metallic 
wall.  An  explanation  of  such  a  process  of  diffusion  may  be 
found  in  Section  D-3  of  this  report. 

It  is  known  that  an  infinite  number  of  modes  of  propagation 
are  found  inside  a  cylinder.  In  order  to  analyze  the  effect  of 
non-uniformity  in  the  most  simple  case,  we  shall  treat  the  case 
where  all  the  components  of  the  electromagnetic  field  are  inde¬ 
pendent:  of  the  coordinate  6,  namely  the  TM^  mode.  The  other 
modes  of  practical  interest  will  be  studied  later. 
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(1)  Analysis  of  the  Problem 

Due  to  the  nature  of  the  propagation  along  the  z-axis 
the  solution  of  the  Maxwell  equations  may  be  written  in  the  form: 
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where  E  and  H  are  functions  of  r.  The  quantity  h  is  the  eigen- 
o  o 

value  which  is  to  be  determined  from  the  boundary  conditions 
which  shall  be  stated  later. 

The  system  of  equations  which  define  the  non-zero  com¬ 
ponents  of  the  electromagnetic  field  is: 
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where  K2  is  defined  by  the  expression 
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According  to  the  discussion  in  Section  D-3  concerning 
the  diffusion  inside  a  cylinder,  the  differential  equation  for 
n^  when  the  plasma  is  radially  non-uniform  is: 


d=ne  !  dne  a 

— j"  2"  +  —  ~ —  +  —  n  =0 
dr  r  dr  D  e 


(D-54) 


This  implies  that  the  electron  density  which  appears  in  the 
Eq.  (D—  53)  can  be  defined,  in  agreement  with  Eq.  ,  v-54)  ,  by  the 
expressions: 


n  =  n  C7  r  if  ^  >  0 

e  eo  o  i  V  D  D 


r  a.  .  n 

n  =  n  I  /—  r  if—  <0 
e  eo  o  v  D  D 


where  n  is  the  electron  density  at  the  axis  of  the  cylinder, 
eo 

The  system  of  Eqs.  (D-50) ,  (D-51) ,  and  (D-52)  can  be 
reduced  to  the  solution  of  the  following  equation: 


i2E  T  Ju2l  d  E 

rA  :-aT  — —  +  (ksK2-h2)E  =  0  (D-55) 

dr2  r  Ks(k2lr-h2)  dr  dr  oz 


In  this  problem  the  z-component  of  the  electric  field  has  to  be 
zero  at  the  wall.  This  condition  will  be  used  to  determine  the 
eigenvalues  of  h  in  Eq.  (D-55) .  Since  it  is  not  possible  to 
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find  a  solution  to  Eq.  (D-55)  in  closed  form,  the  eigenvalues 
must  be  found  by  numerical  procedures . 

(2)  Nondimensional  Analysis 

As  stated  previously,  a  numerical  computation  is  neces 
sary  to  determine  the  eigenvalues  of  h  appearing  in  Eq.  (D-55) . 
In  order  to  do  this,  it  is  convenient  to  transform  Eq.  (D-55)  to 
a  nondimensional  form.  Let  us  define: 


[E  ] 
ozJ 

r  e  ] 

oo  J 


E' 

oz 


[k]  [r]  =  r 


OL 

Dk2 


[  of  1 

[c2] 


[k2] 


Thus,  we  have  used  [k]  in  order  to  reduce  the  lengths  to  a  non- 

dimensional  form  and  E  ,  the  electric  field  at  the  axis  of  the 

oo 

cylinder,  to  reduce  the  Eqz  component  to  a  nondimensional  form. 

If  we  multiply  Eq.  (D-55)  by  ; — — -p-  3-r  ,  we  obtain: 

1  Eoo  i  L  K  J 


d2  E ' 
oz 

dY2 


=  +  g(r) 


d  E1 


~°Z  +  [  K2  -  h2  ]  E '  =0 

dr  oz 


(D- 


56) 
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where 

-h2  .  _*l£ 

g(r)  =  ( K2  -H3  )  K2  Br 


(D-57 ) 


From  these  equations,  the  eigenvalues  of  h  may  be  determined. 
(3)  Determination  of  the  Eigenvalues  of  h 

Equation  (D-53)  can  be  written  in  the  form: 


¥? 


K2 


1  - 


1  - 


ts  xi 

eo 

/_0L_ 

m  e  cJ3  Jo 

,  /  Dk2 

e  o 

e2  n 

f  I - 

eo  ,  j 

/  “ 

m  c  a?  o 

i  v  Dk3 

e  o 


if 


0 


if 


a 

D 


<  0 


The  plasma  frequency  co^2 
formula: 


at  the  axis  is  defined  according  to  the 


e*n 

eo 

m  f 
e  o 


OC  P 

Thus,  there  are  two  parameters:  and  •  These  parameters 
were  used  in  the  numerical  determination  of  h.  k  was  used  as  a 
variable. 


(4)  Results 

In  order  to  present  the  results  in  a  conventional  form, 
CO 


/—  R  and  — ^  are  selected  as  parameters,  where  R  is  the  radius 
VD  o  CO  o 

of  the  cylinder.  The  index  of  refraction  of  the  waveguide. 


defined  according  to  the  following  expression: 
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has  been  plotted  in  Figs.  1-4  as  a  function  of  the  signal  wave 

X  60  2ff 

length  — •  for  four  values  of  — where  X  =  —  . 


60 


It  is  seen  that 


as  X  approaches  zero,  — ~  must  approach  1, 

vf 


In  each  of  the  figures  the  curve  corresponding  to  the 
limit. case  of  w  =  0  is  shown.  This  particular  case  is  defined 
by  the  equation: 


|  2.405\s  Xa 

1  2ir  R  8 
'  o 


1 


I  J 


> 


Fig.  D-4.1 

V  ave  Guide  Index  of  Refraction  as  a  Function  of  Signal  Wave  Length  for 


Wave  Guide  Index  of  Refraction  i  "■  a  Function  of  Signal  Wave  Length  for  =  0.  316 


WiSTBVWf.  N.Y. 
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APPENDIX  I 


SUMMARY  OF  CHEMICAL  REACTIONS  AND  THEIR  RATE  CONSTANTS 


This  appendix  summarizes  a  large  number  of  chemical  reactions 
which  are  significant  at  the  high  temperatures  of  reentry  into 
the  earth's  atmosphere  at  hypersonic  speeds. 

The  presentation  here  includes  two  classes  of  units, 
namely  cm,  sec,  particle  and  meter,  sec,  kg  mol'  (required  for 
the  present  analysis).  The  tables  are  sel  -explanatory,  though 
a  word  of  caution  is  in  order  for  the  ionization  reactions  10 
through  15  The  rate  constants  for  these  reactions  are  tabu¬ 
lated  as  the  forward  (i~e  to  the  right)  values,  which  is 


contrary  to  the  remainder  of  the  appendix. 
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SUMMARY  OF  CHEMICAL  REACTIONS  AND  RATE  CONSTANTS 
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APPENDIX  II 


ONE-DIMENSIONAL  STEADY  STATE  DIFFUSION  ACROSS 
A  DISCONTINUITY  IN  THE  THERMODYNAMIC 
PROPERTIES  IN  AN  IONIZED  GAS* 


For  the  special  case  of  a  one-dimensional  steady-state 
problem  with  the  x-coordinate  (with  respect  10  which  the  non¬ 
uniformity  occurs)  normal  to  a  planar  discontinuity  in  thermo¬ 
dynamic  properties,  the  reduced  electron  and  ion  momentum  equa¬ 
tions  simplify  to  the  forms: 


, .  dn 

d<P  ^  e 
nu  -  n  (1  —  -  D  — — 

3  e  e*e  dx  e  dx 


(II-D 


-  -  dn. 

d<|  ^  i 

n.u.  =  -n.il.  —  -  D,  -r — 
l  i  i  i  dx  l  dx 


(II-2) 


where 


u  =  — ^ —  =  electron  mobility 
e  m  v 
e  e 

4.  =  — ^ -  =  ion  mobility 
i'  i 


kT 

D  =  - -  =  electron  diffusion  coefficient 

e  my 
e  e 


2kT 

D.  =  -  =  ion  diffusion  coefficient 

1  m  .  v  . 

1  1 


*This  Appendix  was  prepared  by  Dr.  Frank  Lane 
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Next,  Eqe.  (II-l)  and  (II-2)  are  divided  by  and  respectively, 
differentiated  with  respect  to  x,  and  added.  This  gives 


_L  A. 

»e  d* 


e  e 


-JL. 

M . 

~  (n. 
dx  i 

u.) 

i 

”  dx 

i 

D 

d2  n . 

D, 

d2n . 

e 

l 

l 

^e 

dx3 

dx2 

n 


d$ 
e  dx 


_d_ 

dx 


d$ 

r'i  dx 


( II-3) 


Next,  the  assumption  that  n_^  is  nearly  equal  to  n^  is  in¬ 
voked  (i.e.  their  difference  is  a  quantity  small  of  first  order). 


Furthermore  the  relation  between  A  $  ana  (n.-u  )  shows  that  the 

1  e 


electric  field 
the  potential 


£x 
d$\ 
dx 


should  also  be  small  of  first  order.  Then 

terms,  whose  difference  occurs  in  (II-3) ,  may 

be  dropped.  Using  the  equivalence  of  the  electron  and  ion  pro- 

ds  ng  d2  n^ 

duction  rates  and  setting  -  ~  dx2"  '  t^iere  results  the  follow¬ 

ing  relation: 


h  I—  +  — 

e  u  u . 

\  e  1 


d2  n  1 


dxc 


D  D.  ' 
e  l 

U  H. 

e  , 


(II-4) 


or 


where 


D  = 


n  = 
e 


D  ji.  +  D.jj 
e  l _ l  e 

u  t-  U . 

e  i 


d2  n  D  a .  +  D,  u 
e  ei  le 


dx2 


H  +  U . 
e  'l 


=  -  D 


d2  n 
_ € 

dxz 


( II-5) 


=  ambipolar  diffusion  coefficient. 
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If  we  now  approximate  the  electron  or  ion  production  rate 
the  expression 

/  n  3  \ 


n 

e 


a 


i 


h  by 
e 


(II-6) 


where  ne  is  equilibrium  electron  number  density  and  a  is  the 
o 

electron  creation  rate  per  unit  volume,  then  Eq.  (II-5)  takes 
the  form 


D 


d2  n 

_ e 

dx2 


+  a 


o 


(II-7) 


This  is  the  equation  which  must  be  solved  subject  to  boundary 
conditions  appropriate  to  the  physical  situation  of  interest. 

In  order  to  get  an  idea  of  the  distances  over  which  elec¬ 
tron  density  is  "smeared  out"  in  making  the  transition  between 
two  asymptotic  values  (such  as  might  be  expected  across  a  shock) 
tne  following  idealized  one-dimensional  situation  is  proposed. 
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A  plane  nterface  is  conceived  across  hich  temperature  jumps 
discontinuously .  The  temperature  is  assumed  constant  in  each 
half  space.  The  electron  density  assumes  its  equilibrium  value 
ne  +  as  x  -»  +  ®,  and  a  different  higher  equilibrium  value  neQ 
as  x  -  ® .  The  system  is  to  be  analyzed  for  electron  density 
distribution  using  Eq.  (II-7)  in  each  half  space  amd  matching 
up  properly  the  two  solutions  at  the  interface  (x  =  0) .  The 
matching  conditions  are  continuity  of  electron  density 


n  (0+)  =  n  (0  ) 
e  e 


(H-8) 


and  continuity  of  electron  current 


n  u 
e  e 


=  n  u 
e  e 


(II-9) 


0 


Now,  from  the  equivalence  of  electron  and  ion  production 
rates,  we  find  that 


-r—  (n .  u  .  -  n  u  )  =  0 
dx  li  e  e 


(11-10) 


This  is  equivalent  to  saying  that  Div (current)  =  0  and  since  cur¬ 
rent  vanishes  at  x  -*  ±  ®  therefore  current  =  0. 


n  u  = 
e  e 


or 


n ,u . 
l  i 


(11-11) 
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Therefore,  utilizing  ( II— 1 )  and  (II-2)  we  cind  that 


n  u 
e  e 


n .  u . 

l  l 


and 

(II-2) 

d$ 

dn 

D  _ § 

dx 

e  dx 

d$ 

dn . 

n  _ i 

dx 

Di  dx 

(11-12) 


(11-13) 


Dividing  (IT-12)  by  /ig  and  (11-13)  by  ,  utilizing  the  equiva¬ 
lence  (11-11)  and  adding,  we  find 


n  u 
e  e 


_L  _L 

+  "i 


\  ^e 


=  (n  -n.)  — 
e  i  dx 


D  D  .  \  dn 
_ e  __i  \  _ e 

"e  +  "i  dX 


(II-14) 


where  again  n^  ~  n^  has  been  invoked  in  the  last  term.  Since 

(n  -n.)  is  small  of  second  order,  (11-14)  leads  to  the  result: 
e  i  dx 


dn 

_  € 

n  u  ~  -  D  — — 
e  e  dx 


(II-15) 


Therefore  the  second  interface  condition  (II-9)  reduces  to 
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dn 

_ c 

dx 


=  D 


dn 
_ _ € 

dx 
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(11-16) 
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Now  the  general  solution  to  the  Eq.  (II-7)  may  be  written  in  the 


form 
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(H-17) 


where  K  is  an  undetermined  constant. 

The  coefficients  oc  and  D  and  the  equilibrium  electron  den¬ 
sity  nQ  (as  well  as  the  temperature)  differ  in  each  half  plane. 
Hence  (11-17)  must  be  applied  individually  to  each  half  plane. 

In  the  right-hand  half  plane  the  (-)  sign  applies  in  the  exponent 
whereas  in  the  left-hand  half  plane  the  (+)  sign  applies.  Also 

the  constant  K  differs  in  the  two  half  planes.  Moreover,  in  the 

cin 

process  of  finding  this  solution  the  following  expression  for  — 
ir  found: 


dn 
_ _ e 

dx 


± 


(11-18) 


where  the  (-)  sign  applies  in  the  right-hand  half  plane  and  the 
(+)  sign  applies  in  the  left-hand  half  plane.  This  now  consti¬ 
tutes  sufficient  information  to  solve  for  the  values  of  K  appro¬ 
priate  to  the  right-hand  and  left-hand  half  planes.  The  actual 
numerical  work  requires  solution  of  a  cubic,  and  this  was  per¬ 


formed  graphically. 
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Following  is  a  table  of  results,  Table  I.  In  the  numerical 

work  n_  r  was  assumed  to  be  106  and  r,_  was  given  the  values 
eo  eo 

107  ,  108 ,  •••  1012.  The  values  in  the  table  were  computed  using 

a  constant  density  p  =  p  atmospheric  x  10-4  and  temperatures 
corresponding  to  this  density  and  the  equilibrium  electron  num¬ 
ber  densities  tabulated.  This,  of  course,  departs  from  jump 
conditions  across  a  shock  but,  nevertheless  the  results  should 
be  indicative  of  the  way  in  which  electron  densities  vary  between 
asymptotic  values. 

/v? 

The  representative  lengths  J  were  computed  by  comput¬ 

ing  a.  and  the  ambipolar  diffusion  coefficient  D  corresponding  to 
the  appropriate  density  and  temperature. 


and 


ess  4  k  T 

D  -  -  since  m.  »  m 

m .  v  ■  i  e 

i  i 


«  2''2 

V.  =  -7=-  S,  n 
1  '  dr 


/  m .  +rr. 

[ _ i_ 

n 

1  m . 

m 

\  1 

n 

(11-19) 


(11-20) 


where 


S „  =  collision  cross  section  for  ions  and 
d 

neutral  particles 


n^  =  neutral  particle  number  density 
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TABLE  I 


1/cc 

■f 

ne 

to 

1/cc 

c 

um 

cm 

i 

NP 

i 

1/ 

v 

T 

+ 

T 

n 

n 

ne0 

v  lot 
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m 

8 

i\ 

+ 

i\ 

(°K) 

(°k: 

n+ 

nn- 

106 

107 

330 

312 

.2272 

.0395 

2175 

2380 

2.48xl016 

2.57xl016 

106 

| 

108 

330 

260 

.5407 

.0638 

2175 

2625 

2.48xl016 

2.68x101b 

1  10° 

109 

330 

100 

.8041 

.0712 

2175 

2990 

2.48xl016 

2.83xl015 

106 

io10 

3  30 

32. 

5 

.9315 

.0724 

2175 

3480 

2.48xl01B 

3 . OOxlO1 5 

106 

1011 

330 

12 

.9774 

.0738 

2175 

4180 

2.48x101b 

3.28xl01B 

106 

1012 

3  30 

6 

.9918 

.0793 

2175 

5270 

2.48x101b 

3.83xl016 

— 2 -  =  10-4 

o 

atraos 

=  neutral  particle  number  density 

m  =  neutral  particle  mass 
n 

and  the  radicand  simplifies  since 


m ,  »  m 
1  n 


Expression  (11-20)  is  derived  from  consideration  of  the  momentum 
transfer  in  a  collision  between  an  ion  arc  a  neutral  particle. 
The  collision  cross  section  is  given  approximately  by 


S,  = 
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r .  +r 
i  n 


2 


\ 


IT 


2 


(H-21) 
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where  a  smooth  sphere  model  is  assumed  and  r.,  r  are  the  ion 

in 

and  neutral  particle  radii  which  are  assumed  to  be  equal  and  of 
the  order  of  1C~S  cm.  The  electron  creation  rate  per  unit  vol¬ 
ume  a.  is  computed  as  follows: 

The  dominant  reaction  is  assumed  to  be 

N  +  0  *  N0+  +  e  (11-22) 


The  rate  of  electron  production  is  given  by 


dn 

— —  =  k 
dt  d 


n 
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n 


k  n .  n 
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(11-23) 


but 
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so  that 


n  =  k,  n 
e  d  n 
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n  n, 
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n  n, 
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k,  n 
d  n  0 


n  2  \ 
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(11-24) 


(H-25) 


where 


n 

n 


n 


C 


k 

r 


equilibrium  electron 
number  density 
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and  where  we  now  see  that 

a.  =  kJ  n  n  (11-26) 

d  n  0 

Now  n^  are  computed  from  tables  (for  air)  at  the  appropriate 
temperature  and  density,  and  k^  is  computed  using  methods  of 
classical  chemical  kinetics. 

Figure  E-2.1  gives  the  actual  computed  electron  density 
distribution  corresponding  to  each  of  the  six  asymptotic  condi¬ 
tions  listed  in  Table  I.  The  transition  becomes  sharper  as  the 
hot-side  electron  density  becomes  higher. 
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APPENDIX  III 

PERIODICALLY  DISTURBED  PLASMA 


INTRODUCTION 

The  simple  criterion  used  to  determine  radar  cross  section 

of  an  ionized  wake  is  to  assume  that  for  electron  density  regions 

W  3  w  3 

where  <  1,  no  reflection  occurs,  and  if  %  "  >  1  total  reflec- 

cf  u) 

tion  occurs.  This  simple  criteria  is  useful  if  the  collision 
frequency  is  low  and  the  density  gradients  are  large  through  the 
above  critical  conditions. 


The  problem  of  interest  here  is  the  electromagnetic  charac- 

0)  3 

teristics  of  *  plasma  region  where  — «  1,  but  whe^e  stationary 

W  o 

a)  2 

or  time  varying  regions  of  higher  density  (  still  less  than  1) 
occur  with  spacing  comparable  to  the  half  wavelength  or  multiple 
of  the  incident  radiation.  This  situation  is  comparable  with  a 


multilayered  dielectric  filter,  except  that  in  general  the  loca¬ 
tions  of  the  maximum  electron  density  regions  are  randomly  changing. 

It  must  be  noted  that  while  it  is  possible  to  obtain  a  strong 
reflection  from  such  a  periodic  dielectric  structure,  it  is 
equally  probable  that  total  absorption  will  occur,  by  the  same 
process.  However,  if  the  stxu.cr.ure  is  randomly  varying,  the 
process  of  interest  is  that  of  reflection,  which  ultimately  is 
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detected.  The  totally  absorbed  energy  is  equivalent  to  a  trans¬ 
parent  region  as  far  as  the  reflected  characteristics  are 
concerned . 

EXPERIMENTAL  OBSTACLE 

In  order  to  properly  formulate  the  analytical  model  required 

to  generalize  the  problem,  some  basic  experimental  information 

must  be  supplied.  Of  course  such  information  must  be  sufficiently 

quantitative  that  a  true  model  can  be  formulated.  Therein  lies 

the  experimental  problem,  for  since  the  basic  plasma  required 

03  2 

is  strongly  underdense,  — ^  «  1,  measurement  of  the  plasma 

03 

properties  becomes  extremely  difficult. 

One  approach  to  this  measurement  problem  is  to  induce  the 
required  perturbation  in  such  a  manner  that  the  plasma  properties 
can  be  accurately  deduced  from  other  parameters. 

One  such  experimental  arrangement  considered  was  to  propa¬ 
gate  a  strong  sound  wave  into  the  shocked  gas.  With  such  an 
arrangement-  the  amplitude  and  spacing  of  the  periodic  structure 
could  easily  be  changed  within  limits.  Since  the  maximum  ex¬ 
pected  variations  in  pressure  would  be  small  compared  to  the 
base  pressure  the  electron  density  variation  might  be  assumed 
to  be  directly  proportional  to  the  pressure  variation.  That  is 
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While  all  parts  of  this  scheme  are  correct  in  concept, 
severe  problems  arise  in  coupling  the  required  sound  energy  to 
the  shock  tube.  Although  several  different  coupling  systems 
have  been  tried,  the  attempt  to  use  sound  waves  to  induce  plasma 
variations  was  abandoned  because  of  the  extensive  facility 
+-hat  would  be  required. 

To  obtain  an  order  of  magnitude  estimate  of  the  effect  of 
a  turbulent  wake  on  the  reflected  radiation,  an  experiment  was 
performed  in  the  shock  tube.  A  cylindrical  bar  was  mounted 
across  the  tube  parallel  to  the  direction  of  polarization  of 
the  radiation.  Several  diameters  were  tried  up  to  V*  which 
was  the  maximum  diameter  consistent  with  an  unchoked  flow. 
Underdense  plasma  conditions  were  run  and  observations  were 
made  of  the  reflected  signal  with  the  use  of  the  electric  field 
detector,  in  the  polarized  and  antipolarized  directions. 

Even  though  the  scale  of  the  turbulence  expected  was  an 
order  of  magnitude  smaller  than  che  wavelength,  if  vex'y  strong 
effects  were  to  be  associated  with  this  type  of  wake,  a  noise 
signal  would  be  imposed  on  the  reflected  power  signal.  Since 
no  such  signals  were  observed,  it  may  be  concluded  that  small 
scale  effects  are  unimportant  in  this  process. 
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RESULTS 

An  attempt  was  made  to  see  if  strong  effects  could  be  ob¬ 
tained  under  conditions  where  the  scale  of  the  disturbance  was 
comparable  to  the  wavelength,  even  if  no  precise  description  of 
this  disturbance  could  be  given. 

The  basic  scheme  used  was  the  addition  of  a  "Tee"  section 
to  the  shock  tube,  where  the  area  of  the  arm  of  the  "Tee"  is 
about  half  the  area  of  the  tube.  As  the  flow  passes  the  "Tee" 
section,  under  certain  conditions,  a  "whistle"  type  effect  can 
occur.  If  this  is  the  case,  periodic  pressure  fluctuations 
should  be  induced  in  the  flow.  Such  pressure  variations  should, 
in  turn,  create  regions  of  higher  electron  density. 

In  Fig.  E-3.1  a  trace  is  presented  showing  the  pressure 
fluctuations  occurring  in  the  arm  of  the  "Tee"  near  the  main 
flow.  Although  such  variations  could  only  be  obtained  for  a 
very  restricted  set  of  operating  conditions.  Fig.  E-3.1  clearly 
shows  that  some  type  of  periodic  pressure  producing  mechanism  is 
occurring. 

To  check  the  validity  of  the  result,  a  simple  step  model 
for  the  plasma  behind  the  shock  can  be  used  to  compute  the  re¬ 
flected  signal.  Assuming  the  pressure  defect  is  carried  with 


the  flow,  the  simple  model  for  the  plasma  postulates  a 
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reflection  from  the  front  of  the  shock  and  a  second  reflection 
from  the  pressure  disturbed  region  which  exists  behind  the  shock 
a  distance  (z)  where  z  =  v  t  with  the  particle  velocity  with 
respect  to  the  shock  and  time  t  starting  from  the  creation  of 
the  disturbance. 

For  the  moment,  if  it  is  assumed  that  the  waveguide  wave¬ 
length  is  not  appreciably  altered  in  the  first  plasma  region 

directly  behind  the  shock,  then  it  is  clear  that  when  the  dis- 

X 

tance  z  the  reflected  signal  from  the  disturbed  region  will 

be  180°  out  of  phase  from  the  reflection  occurring  at  the  front 
of  the  shock,  thus  producing  a  minimum  in  the  standing  wave  ratio. 
On  the  other  hand,  for  z  =  X^/2  the  two  reflected  signals  are  in 
phase  and  a  maximum  should  occur  in  the  SWR. 

In  Fig.  E-3.2  a  trace  of  the  electric  field  is  shown  for  the 
critical  whistle  conditions  described.  While  this  result  must  be 
considered  tentative,  it  is  clear  that  before  the  shock  enters 
the  "Tee"  section,  the  SWR  is  small  indicating  an  "underdense" 
plasma.  After  the  transition  of  the  "Tee"  where  the  effects 
cannot  be  clearly  described  at  the  present  time,  the  amplitude 
modulation  of  the  SWR  is  seen  to  decrease  and  then  subsequently 
to  increase.  The  end  transition  occurs  as  first  the  plasma 


passes  the  detector  and  then  reaches  the  antenna,  where  the 
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signal  drops  to  zero.  Since  this  type  of  modulating  behavior  is 
consistent  with  the  simple  model  and  inconsistent  with  other  pos¬ 
sible  mechanisms,  it  is  felt  that  this  trace  represents  the  type 
of  strong  reflections  that  can  be  obtained. 
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Top  Trace:  Heat  transfer  gauge 
located  at  entrance  to  'Tee" 
section. 

Lower  Trace:  Pressure  in  arm 
of  "Tee". 

Sweep  Time:  100  (is/cm. 


Figure  E  -3.  1 


Electric  Field  Detector.  Located 


3  meters  downstream  of  "Tee' 


section. 


Sweep  Time:  200  jis/cm. 


Note:  The  cancellation  and 


subsequent  growth  of  the  signal  as 
the  "filter”  structure  goes  through 

Xg 

the  — “■  condition. 
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